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ABSTRACT 


This  document  provides  the  final  report  for  Contract  Number 
N00167-82-C-0023,  "Applications  of  an  Improved  Wave  Drag  Prediction  Method." 
The  method  is  based  on  a  modified  supersonic  linear  theory  computer  technique 
developed  in  this  contract  and  two  earlier  contracts,  N00167-78-C-0005  and 
N00167-79-C-0123,  all  of  which  were  jointly  funded  by  DTNSRDC  and  NASA/Langley 
Research  Center.  The  modified  linear  theory  differs  from  ordinary  linear 
theory  in  that  it  uses  the  exact  boundary  condition,  it  uses  the  local 
perturbed  velocity  to  calculate  s  =/m^-1 ,  it  uses  the  exact  pressure 
coeffficient  equation,  and  it  uses  characteristic  tracing  to  determine  regions 
of  influence.  The  theoretical  and  numerical  techniques  are  described  in 
detail.  A  computer  user's  manual  is  included  which  provides  a  detailed 
description  of  the  inputs  and  outputs  of  the  computer  routine,  including 
recommended  guidelines  for  preparing  the  geometric  input.  Comparisons  between 
predictions  of  drag  or  pressures  using  the  modified  linear  theory  are  made  for 
five  cases  with  wind  tunnel  data  or  other  theories.  Tnese  comparisons  verify 
the  improved  accuracy  of  the  modified  linear  theory  relative  to  ordinary 
linear  theory,  particularly  for  high  rlach  numbers  and  high  surface  slopes. 


IMPROVED  WAVE  DRAG  PREDICTIONS  USING  MODIFIED  LINEAR  THEORY 

1.0  INTRODUCTION 


Supersonic  linearized  theory,  including  the  special  cases  of  slender  body 
theory  and  area  rule,  have  well  known  limitations.  Uhile  linear  theory  does 
remarkably  well  in  predicting  lift  and  pitching  moment,  it  does  not  do  as  well 
in  predicting  drag.  The  largest  discrepancies  occur  when  the  sweep  angle  is 
nearly  equal  to  the  Mach  angle  (cos_ll/M).  Inaccurate  prediction  results 
because  drag  is  the  integral  of  pressure  times  local  slope,  and  the  largest 
errors  in  predicted  pressure  are  likely  to  occur  where  the  slopes  are 
largest.  Van  Dyke's  second-order  method  gives  much  improved  accuracy  relative 
to  slender-body  or  linear  theory  results  for  two-dimensional  flow,  either 
planar  or  axisymmetric,  but  a  method  for  extending  it  to  three  dimensions  has 
not  been  found.  For  aircraft  configurations,  analytical  methods  for 
predicting  drag  more  accurately  than  area  rule  have  not  been  found,  short  of 
the  very  complex  "exact"  solutions  such  as  method  of  characteristics, 
time-dependent  equations  of  motion,  or  relaxation  techniques.  These  exact 
methods  require  large  amounts  of  computer  time  and  are  often  sensitive  to 
input  data  smoothness,  choice  of  arbitrary  parameters,  etc.  These  techniques 
are  not  ideal  for  preliminary  design,  where  many  answers  are  required 
quickly.  Preliminary  design  requires  a  method  which  approximates  linear 
theory  in  complexity  and  approaches  the  exact  solutions  in  accuracy.  The 
modified  linear  theory  technique  has  been  developed  to  provide  an  efficient, 
accurate  method. 

The  analytical  and  numerical  techniques  used  in  the  modified  linear  theory 
analysis  routine  are  described.  Calculations  for  several  cases  are  compared 
•with  wind  tunnel  data  or  other  theoretical  methods.  A  computer  user's  manual 
is  included  as  Appendix  B. 

2.0  LINEAR  THEORY  PERTURBATION  RELATIONS 


Because  line  distributions  of  sources  and  sinks  are  inadequate  for  a 
general  three-dimensional  solution,  all  calculations  described  in  this  paper 
utilize  surface  distributions  of  sources  and  sinks.  The  velocity  potential 
equation  is 
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where  S  is  the  souce  strength,  ds  =  / on  the  surface  of  the 

body,  and  the  integral  is  taken  over  that  portion  of  the  surface  included  in 
the  Mach  forecone  from  the  point  x,y,z.  Now,  under  certain  nonrestrictive 
conditions,  the  perturbation  velocity  components  can  be  derived  to  be  of  the 
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where  e  is  the  local  slope,  vn  is  the  perturbation  component  perpendicular 
to  the  freestream  and  normal  to  v«  and  vp  is  the  perturbation  component 
perpendicular  to  the  freestream  ana  parallel  to  the  local  surface.  The 
functions  F,  G,  and  H  are  dependent  on  the  value  of  b  ,  the  limits  of 
integration,  the  geometry  of  the  model,  and  the  functions  used  to  describe  the 
variations  of  source  strength  in  the  y,  z  directions. 

3.0  MODIFICATIONS  TO  LINEAR  THEORY 

The  accuracy  of  the  current  method  results  from  the  combined  effects  of 
the  exact  boundary  condition,  and  use  of  the  local  (perturbed)  value  of 
6  =■  y,  -  -  ;  .  The  exact  pressure  coefficient  equation  is  also  used. 

The  exact  boundary  condition  (Fig.  1)  requires  not  only  the  use  of  the 
perturbed  streamwise  velocity  component  but  also  the  determination  of  the 
surface  slope  (ei  in  the  plane  defined  by  the  freestream  velocity  vector  and 
the  normal  to  local  surface.  This  three-dimensional  determination  of  the 
slope  and  the  velocity  component  boundary  conditions  is  necessary  even  in 
quasiplanar  cases,  such  as  for  wings. 

The  primary  effect  of  the  local  e  (Fig.  2}  is  on  tne  constant  of 
proportionality,  1/s.  The  1 /e  factor  appears  in  every  term  but  one  on  the 
riynt  side  of  Eqs.  (2-4).  Thus,  it  has  a  direct  effect  on  each  of  the 
perturbation  velocities.  As  the  local  .'lach  number  approaches  l.u,  the  value 
of  1/8  approaches  infinity.  This  is  obviously  an  undesirable  result. 
Therefore,  a  correlated  local  Mach  or  b  has  been  developed  based  on 
calculations  for  two-dimensional  ramps  and  cones.  Tne  correlation  puts  a 
limit  on  the  maximum  value  of  1/b  and  causes  the  correlated  value  to  approach 
the  limit  more  slowly  than  with  the  exact  equation.  The  effect  of  local  3  on 
the  region  of  influence  (characteri Stic  directions)  is  less  pronounced, 
particularly  for  two-dimensional  cases,  either  planar  or  axisymmetric. 

However,  for  three-dimensional  flow  the  variations  in  propagation  direction 
can  be  important  because  of  the  lateral  or  longitudinal  displacement  of 
interference  effects.  Region  of  influence  effects  will  be  discussed  further 
below. 

The  two  modifications  have  opposite  effects.  For  compressions,  the  exact 
tangency  condition  tends  to  reduce  the  magnitude  of  the  perturbations  because 
(1  +  <?■,  /V0  )  is  less  tnan  unity,  and  the  required  normal  velocity  vn  is 

less  than  with  the  linearized  boundary  condition.  But  the  local  b  in  Eqs. 
(2-4)  tends  to  increase  the  perturbations  values  because,  when  M  is  smaller 
than  the  freestream  value,  1/b  is  larger.  Conversely,  for  expansions,  (1+^, 
/V o)  >  1,  and  the  required  vn  is  larger  than  the  linearized  values,  but 
1/b  is  smaller  than  the  freestream  value. 

Theoretical  pressures  predicted  by  the  modified  method  are  compared  with 
exact  theory  and  linearized  theory  for  45-deg  swept  ramps  (Fig.  3).  The 
modified  linear  theory  shows  greatly  improved  correlation  with  exact  theory 
relative  to  ordinary  linear  theory  for  both  expansions  and  compressions. 

Also,  where  the  lineary  theory  predicts  infinite  perturbations  at  sonic  edge 
conditions  for  the  swept  ramp  in  Fig.  3  (M0  =  Jz)  ,  tne  modified  theory  on 
the  expansion  side  shows  no  tendency  toward  infinity.  Use  of  tne  local  3 
eliminates  expansion  singularities  because  1/b  approaches  zero  with  finite 
positive  perturbations.  On  the  compression  side,  the  exact  and  modified 
theories  cannot  be  computed  at  M0  =  /~Z  for  this  two-dimensional  case 


2 


because  the  flow  could  not  actually  remain  two-dimensional.  However,  it  can 
be  seen  that  the  modified  theory  will  not  predict  infinite  perturbations  here 
either,  because  as  <7*., /V„ — *■  -1,  exact  boundary  conditions  cause  vn  and  the 
source  strength  to  approach  zero.  Thus,  an  equilibrium  must  be  reached  with 
the  value  of  <7^/^  between  0  and  -1.0. 

4.0  ANALYTICAL  ANO  NUMERICAL  TECHNIQUES 

An  important  numerical  technique  which  allows  the  efficient  evaluation  of 
the  perturbation  velocity  integrals  in  Eqs.  ( 2 >  -  (4)  is  the  use  of  Simpson's 
Rule  combined  with  stepping  off  from  every  integration  limit  by  0.07  times  the 
step  size.  For  example,  the  exact  value  of  the  integral 

O 

/ 


i  s  n  =  3.1416.  Even  though  there  is  a  square  root  singularity  at  both  limits, 
a  two  step  (minimim)  Simpson's  Rule  integration  with  the  0.07  offset  gives 


./.C't  -'J) 


I 


which  is  accurate  to  0.2  .‘L  The  0.07  offset  is  used  at  all  times.  Thus,  no 
special  treatment  is  requireu  for  singularities,  and  the  offset  is  small 
enough  that  its  effect  on  non-singular  integrands  is  negligible.  This 
technique  was  derived  by  the  author  for  the  square  root  and  logarithmic 
singularities  which  frequently  occur  in  linearized  aerodynamic  theory. 

A  brief  discussion  of  the  reasons  for  choosing  the  format  of  Eqs.  (1)  - 
(4i  is  in  order.  First,  source  distributions  are  used  rather  than  higher 
order  si ngulari ties  because  they  are  simpler,  and  because  they  are  adequate  to 
describe  the  flow  around  bodies  with  thickness  except  possibly  at  thin 
trailing,  side  or  subsonic  leading  edges  where  the  upper  and  lower  surface 
pressures  are  not  equal.  The  assumed  source  strengths  vary  linearly  in  both 
directions  on  a  panel  and  are  continuous  across  panel  boundaries  unless  the 
slope  is  discontinuous.  This  continuity  of  source  strengths  and  the  fact  that 
the  panels  can  be  curved  in  both  directions  allows  the  use  of  larger  (i.e., 
fewer)  panels  to  describe  a  configuration  than  would  otherwise  be  possible. 

Second,  the  form  of  Eqs.  (2)  -  (4)  was  chosen  for  the  following  reasons. 
The  velocity  components  must  be  obtained  from  derivatives  of  the  potential 
function.  The  double  integral  in  the  potential  equation  cannot  readily  be 
evaluated  analytically  and  at  least  one  integration  must  be  performed 
numerically.  Numerically  evaluated  derivatives  are  highly  unreliable, 
particularly  for  the  derivative  of  a  numerically  evaluated  integral.  A  much 
more  reliable,  accurate  and  efficient  method  is  to  perform  one  integration 
analytical ly ,  differentiate  analytically  for  the  three  velocity  components, 
and  then  perform  the  remaining  integration  numerically  for  each  velocity 
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component.  Having  selected  this  method,  the  question  remains  of  which 
integration  to  perform  analytically,  the  streamwise  or  lateral?  both  options 
were  examined.  It  was  found  that  analytically  integrating  laterally  gave  a 
simpler  integrand  with  fewer  singularities  and  was  much  more  compatible  with 
the  translations  required  for  characteristic  tracing.  The  derivations  of  Eqs. 
(2)  -  (4)  for  laterally  curved  and  flat  panels  are  in  Appendix  a. 

Because  the  perturbation  velocities  are  strong  functions  of  the  local 
6-  Jm  2  -  /  ,  and  the  local  b  is  a  function  of  the  perturbed 

velocity,  an  iteration  is  required.  For  the  first  pass,  a  rough  estimate  of 
the  perturbation  velocities  and  b  is  made,  based  on  the  local  slopes.  This 
provides  faster  convergence  than  using  freestream  b  as  the  first  guess,  a 

marching  solution  scheme  is  used.  The  iteration  is  performed  at  the  front  of 

the  configuration,  then  at  the  mid-points  of  the  first  row  and  then  at  the  aft 
end  of  the  first  row.  At  the  front  of  succeeding  rows,  values  of  so  e 
strength  and  s  from  the  aft  end  of  the  first  row  are  used  with  intei  .ation 
as  needed  if  the  slopes  are  continuous.  If  slopes  are  not  continuoi  across 
the  row,  two-dimensional  or  conical  solutions  are  iterated  upon.  Th  ‘he 

solution  is  stepped  back  to  the  mid-points  of  each  row,  etc.  This  p  jss 

makes  the  iterations  more  stable,  as  points  on  the  back  rows  are  no*  '■Curbed 
by  variations  in  local  b  and  source  strength  on  preceding  rows. 

Variations  in  panelling  techniques  were  investigated  and  two  significant 
effects  were  found:  lateral  curvature  and  continuity  between  panels-  The 
panel  edges  must  be  defined  identically  from  one  panel  to  the  next  to  avoid 
differences  in  and/or  negative  square  roots  in  the  hyperbolic  radius 


v  ,  r  .  '  ' 

calculations.  Differences  in  the  hyperbolic  radius  values,  particularly  if 
near  the  singularity  or  zero  value,  can  cause  large  errors  in  the  integrand 
because  the  net  effect  of  all  the  panels  is  a  small  differences  of  large 
numbers.  The  following  conclusions  regarding  curvature  effects  are  based  on 
the  format  of  Eqs.  (2)  -  (4)  and  might  be  different  if,  for  example,  the 
streamwise  integration  were  performed  analytically.  Longitudinal  curvature  is 
not  important;  i.e.,  if  the  exact  value  of  slope  is  used  in  the  boundary 
condition  and  no  gaps  are  artificially  introduced,  the  exact  ordinate  or 
streamwise  curvature  has  a  negligible  effect.  Lateral  curvature  is  important, 
however.  In  the  vicinity  of  the  receiving  point  and  along  the 
characteristics,  it  is  necessary  to  use  the  curved  panel  solution  from 
Appendix  A.  Each  panel  is  approximated  by  a  circular  arc  segment.  The 
curvature  effect  is  used  for  at  least  two  full  panel  widths  away  from  each 
characteristic.  However,  when  x^  <_  X*,  the  flat  panel  analysis  is  entirely 
adequate.  X*  is  defined  as  the  streamwise  location  where  the  complete  circle 
corresponding  to  the  circular  arc  approximation  to  a  given  panel  lies  entirely 
within  the  forecone  from  the  receiving  point. 

The  curvature  effect  is  used  only  for  the  laterally  constant  source 
strength.  The  curved  panel  solution  with  lateral  source  strength  variation, 
besides  being  quite  complex,  gave  undesirably  large  contributions  in  some 
cases.  The  flat  panel  solution  for  lateral  source  strength  variation  appears 
to  be  more  accurate  and  reliable,  and  so  it  is  used  even  on  curved  panels. 

Two  values  of  s  are  used  in  the  method.  The  true  local  b  is  used  to  trace 
the  characteristics,  or  region  of  influence  boundaries.  The  true  local  liach 
and  b  are  computed  from  the  exact  isentropic  relation: 
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The  correlated  s  is  used  in  the  perturbation  velocity  equations.  The 
correlated  b  is  determined  by  solving  for  the  value  which  would  make  the 
modified  linear  theory  agree  with  a  two-dimensional  Prandtl -Meyer  corner 
solution  with  the  same  (V/V0).  From  the  modified  linear  theory  for  a 
two-dimensional  corner, 

t  *  --F- 

f  *  a  (t  - 
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Now  the  Prandtl -Meyer  angle  is 
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and  for  y  =  1.4, 
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or 
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The  exact  Prandtl -Meyer  corner  solution  is: 


v  . 


-aV^  =  _  dx  =  V  sin  (v-vg) 

4VX  =  dx  -  V  COS  (v-Vg)  -  Vg 

Substituting  Eq.  (13)  into  (9),  we  get 

V  sin  (v-v0  =  s  [v  cos  (v-vg)  -  V0] 

or  3corr.  =  sin  (wo) 

COS  (v-vO)  -  Vg/V 


where  both  the  numerator  and  denominator  are  positive  for  expansions  (V  >  Vo) 
and  both  are  negative  for  compressions  (V  <  Vo).  The  Prandtl -Meyer  solution 
is  exact  for  expansions,  and  it  is  also  the  exact  isentropic  solution  for 
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compressions.  The  difference  between  pressure  coefficients  on  two-dimensional 
ramps  with  and  without  shock  effects  are  less  than  2%for  values  of  Cp  <  0.5 
and  M0  <  2.6  as  shown  in  Figure  4.  Eq.  (15)  is  used  to  evaluate  the 
correlated  b  except  wnen  8corr.  >10  it  is  set  to  10,  and  when  true  local 
Mach  number  is  less  than  1.4,  a  smooth  transition  is  made  from  Eq.  (lb)  at 
‘•'X  =  1.4  to  3corr  =  0.5  at  My  <  l.u. 


Evaluation  of  test  cases  on  swept  ramps  indicate  that  much  better  accuracy 
is  obtained  when  the  characteristics  are  not  yawed  even  though  there  is  a 
lateral  velocity  component.  But  the  1 ateral  velocity  component  should  be 
included  in  evaluating  V/V0  for  the  s  calculations. 

It  is  not  clear  whether  the  s  used  in  the  integrands  should  be  the  value 
at  the  receiving  point,  at  the  sending  points,  or  some  average.  In  most 
cases,  the  difference  was  found  to  be  negligible  so  the  receiving  point  value 
was  used.  But  finally,  a  case  was  found  where  the  choice  of  b  made  a 
difference.  On  a  cockpit-like  forebody,  where  the  slopes  were  approaching 
zero  after  being  negative,  the  method  predicted  pressure  much  too  high.  This 
was  corrected  by  using  an  averaged  b  in  the  integrands  defined  by 

-I - 1 _  j  , 

j  ~ 


where  Xs  is  the  X  of  the  receiving  point,  Xs  is  the  X  of  the  sending 
point,  and  e:_  and  3,<  are  the  local  correlated  s's  on  the  left-running  and 
right-running  characteristics,  respectively,  dse  of  the  averaged  a  eliminated 
the  overprediction  of  pressure  behind  the  canopy  and  had  no  effect  elsewhere. 

The  relative  advantages  of  characteristic  tracing  versus  the  usual  linear 
theory  see-through  integration  boundaries  will  be  examined.  A  typical  plot  of 
the  integrand  in  Eq.  (2)  for  an  axi syinmetri c  body  is  shown  in  Fig.  5.  Tne 
singularity  at  X'  is  first  order  (i.e.,  of  the  form  1/xj  and  thus  produces  an 
infinite  result  for  integration  from  one  side  to  the  singularity,  but  the 
i ntegrand  is  equal  and  opposite  at  X’x  =  x*  +  6  and  =  x*  -  <5. 

Therefore,  tne  result  of  the  total  integration  is  finite  if  the  source 
strength  is  continuous.  But,  if  the  body  has  a  slope  discontinuity,  the 
source  strength  is  discontinuous,  and  at  a  point  located  behind  the  corner 
such  that  x*  is  equal  to  the  x  of  the  corner,  the  perturbation  velocities  are 
theoretically  infinite.  It  is  obviously  incorrect  to  get  an  infinite 
perturbation  from  an  area  which  actually  cannot  influence  the  point  under 
consideration.  If  characteristic  tracing  is  performed,  the  point  at  x^  =  x* 
and  u  =  i  will  be  outside  of  the  influencing  region.  With  the  integration 
taken  only  to  the  characteri stic  line,  the  integrand  is  finite  everywhere. 

The  equations  for  the  rate  of  lateral  translation  of  the  characteristics 
or  the  forecone  are  presented.  Let  y/  be  the  local  surface  lateral  direction 
at  tne  characteri stic ,  and  z'  the  local  surface  outward  normal  direction.  For 
unyawed  characteristic  tracing,  we  get  the  simple  result: 

dy,  =  *_  J\  + 

dx,  3 
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For  the  see-through  case,  the  rate  is: 


for  X  > 


for  X  =  Xi 

where  y  and  l!  are  the  coordinates  of  the  receiving  point  relative  to  the 
sending  panel.  The  see-through  rate  is  derived  from: 

I  '  '  V‘_  '  <  ■  ’*  -  o  —  r  £  .,T  "  r  Y.  -  O 

■!  t  /  ■  .a 

The  usual  linear  theory  (or  see-through)  integration  boundaries  are 
compared  with  the  actual  characteristic  boundaries  in  Figure  6.  For  point  A, 
the  difference  in  the  influencing  region  is  not  great.  At  point  B,  the 
see-through  method  is  obviously  unusable  with  surface  sources  because  the  Mach 
forecone  does  not  include  the  surface  anywhere  near  point  8.  This  condition 
occurs  'when  8  e  <_  -  1,  and  while  extreme,  it  does  occur  on  a  fineness  ratio  4 
circular  arc  body  at  Hach  numbers  of  2.0  and  above.  The  method  is  also  usable 
when  8  e  >_  +  1. 

As  shown  in  Reference  1,  the  characteristic  tracing  method  is  more 
accurate  for  moderate-to-hi gh-angl e  cones,  while  the  see-through  method  is 
better  at  predicting  the  maximum  expansion  pressures  on  slender  circular  arc 
bodies. 

The  characteri stic  tracing  technique  eliminates  contri butions  from  areas 
which  cannot  influence  the  receiving  point,  such  as  a  strut  or  pod  on  the 
lower  side  of  a  wing  which  cannot  affect  the  upper  surface  of  tne  wing. 
However,  for  multi-bodies  and  bodies  with  step  discontinuities  such  as  inlets, 
the  characteristic  tracing  adds  considerable  complexity  to  the  computer  logic 
required.  The  complexity  arises  because  of  multiple  possibilities  for 
determining  regions  of  influence.  Vihen  a  receiving  point  is  downstream  of  a 
step,  the  region  of  influence  of  panels  upstream  of  the  step  can  be  determined 
by:  (a)  tne  tracing  through  space  of  a  characteri stic  from  anywhere  along  the 
step,  (b)  the  tracing  of  a  characteristic  from  the  edge  of  the  step  along  the 
surface,  or  (c)  the  tracing  of  a  characteristic  along  the  surface  from  any 
location  defined  by  (a)  above.  In  addition,  these  complex  regions  of 
influence  can  be  different  for  different  receiving  points,  so  the  calculation 
can  not  be  done  once  and  used  repeatedly. 

It  is  concluded  that  the  characteristic  tracing  method  is  superior  to  the 
see-through  because  of;  1)  the  elimination  of  the  infinite  perturbation  from  a 
slope  discontinuity;  2i  the  capability  to  calculate  accurately  when  s  t  <_  -  1 
and  to  calculate  with  e  e  >  1;  3)  better  numerical  accuracy  on  moderate-and 

high-angle  cones;  and  4)  tTie  elimination  of  contributions  from  areas  which 
cannot  actually  influence  the  point  being  considered.  These  benefits  outweigh 
the  overprediction  of  maximum  expansion  pressures  on  slender  circular  arc 
bodies  and  the  additional  complexity  of  the  calculations. 

Having  selected  the  characteristic  tracing  method,  translation  of  the 
receiving  point  coordinates  is  required  to  maintain  proper  values  of  the 
hyperbolic  radius.  For  axisymmetric  bodies,  the  only  translation  required  is 
an  x  translation.  For  correct  results,  it  was  found  necessary  to  keep  the 
hyperbolic  radius  value  equal  to  zero  on  the  characteri sties  in  the  vicinity 


dy/  =  X-XL  . 

dx,  ~  n- 

y 

dy,  _  =  +  A  -  &2C2 

dx7  b 
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of  the  receiving  point.  By  numerical  experimentation,  this  "vicinity"  was 
defined  by  the  angular  difference  between  the  y  -  z  projection  of  the  surface 
normal  at  the  receiving  point  and  the  normal  on  the  characteristic  being  less 
than  cos  (0.16)  =  80.8  .  After  a  characteri stic  has  left  the  “vicinity"  of 
the  receiving  point,  the  x  margin  from  the  forecone  (hyperbolic  radius  =  0)  to 
the  characteristic  is  allowed  to  become  positive  but  it  is  not  allowed  to 
decrease.  I.e. ,  As  the  characteristics  are  traced  forward,  x  translation  will 
be  done  if  necessary  so  that  the  margin  is  always  greater  than  or  equal  to  the 
margin  at  the  previous  point  on  the  same  characteristic.  The  translated 
receiving  point  coordinates  (Qx,  Qy ,  Qz)  are  single-valued  functions  of  xj, 
the  x  integration  variable. 

For  non-axi symmetric  bodies,  y  and  z  translation  may  be  required.  The  y, 
z  shifting  will  be  considered  in  two  parts.  The  lateral  part  is  in  the 
direction  parallel  to  a  line  between  the  two  characteristics,  and  the  vertical 
part  is  perpendicular  to  the  same  line. 

The  lateral  translation  is  required  on  general  body  shapes  so  that  the 
hyperbolic  radius  can  be  zero  on  both  characteristics  in  the  vicinity  of  the 
receiving  point,  and  to  maintain  the  proper  margins  or  values  when  one  or  both 
characteristics  are  not  in  the  vicinity. 

Vertical  translation  is  required  to  prevent  bulges  such  as  canopies  from 
projecting  into  or  near  to  the  forecone  from  a  point  behind  the  canopy.  This 
is  not  a  problem  caused  by  characteristic  tracing;  rather,  it  can  be  handled 
because  of  characteri stic  tracing.  At  a  given  xi,  all  panel  edges  between 
the  characteristics  are  tested  and  Qy  and  Qz  are  shifted  according  to  the  most 
critical  panel  edge  and  the  criteria  z C£  1  rl/3. 


where  r^  is  the  radius  of  the  circle  passing  through  the  panel  edge  and  the 
two  characteristics,  and  r2  is  the  radius  of  the  forecone  at  X^. 

On  some  general  body  shapes,  the  use  of  characteristic  tracing  results  in 
an  anomaly.  This  anomaly  cannot  occur  on  axi symmetric  bodies,  but  will  occur 
on  shapes  such  as  elliptic  cones  or  bodies  with  flat  or  concave  surfaces,  for 
example.  When  the  inward  normal  from  a  characteristic  trace  has  passed 
through  the  y-z  projection  of  the  receiving  point,  the  value  of  hyperbolic 
radius  increases  as  the  characteristic  is  approached  laterally  instead  of 
decreasing  or  approaching  zero.  Usually,  on  a  panel  cut  by  a  characteristic, 
the  integrand  has  the  opposite  sign  relative  to  neighboring  panels.  But  when 
the  hyperbolic  radius  increases  to  the  characteristic,  the  integrand  has  the 
same  sign  as  neighboring  panels  and  larger  magnitude,  resulting  in  a  large, 
erroneus  contribution  from  panels  which  are  on  the  opposite  side  of  the  body 
from  the  receiving  point.  This  large  undesirable  contribution  has  been 
eliminated  by  neglecting  the  integrand  when  the  above-mentioned  condition 
exists  and  a  lateral  limit  is  a  characteristic. 
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The  non-integral  terms  in  Eqs.  (A-21),  (A-22),  (A-29)  and  (A- 30)  occur  at 
the  receiving  point  when  the  point  is  on  the  body,  and  they  also  occur  in  most 
cases  when  the  forecone  intersects  a  nearby  body  or  when  two  panels  Intersect 
laterally  at  approximately  a  right  angle. 


Consider  a  single  panel  with  a  side  edge  at  ys. 


yfT 
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When  y>yft,  the  non-integral  term  is  a  unit  value  and  is  the  entire  effect, 
because  the  integrands  are  zero.  When  y>ys,  the  non-integral  term  is  also 
the  unit  value,  but  at  y=ys  it  is  half  of  the  unit  value,  and  when  y<y$, 
it  is  zero.  The  integral  contribution  is  such  that  the  net  effect  vanes 
smoothly  from  the  unit  value  at  y=y^  through  one-half  at  y=ys  to  zero  at 

y=ye. 


non-integral  term 


!  •> 

,'r''  ’'—net  effect 
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It  was  found  that  the  numerical  integrals  in  Eqs.  (A-29)  and  (A-3U)  could  be 
significantly  in  error  at  a  lateral  corner  because  the  term  (X-Xj)/^1  - 

(Z-Z^)2  in  the  denominator  becomes  very  small.  In  the  vicinity  of  the 

characteristic,  let  \  *  \  -  t,  where  ( X-XS^%MZ-ZS)2  - 

(y-ys)2  =  0.  Then  the  singular  and  near-singular  ten, is  are  of  the  fonn: 
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If  we  neglect  the  second  order  terms  in  t,  we  have: 
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If  we  apply  a  ratio  to  the  Simpson's  Rule  contribution  at  the  usual  .U7 
offset,  and  set  the  Simpson's  Rule  result  equal  to  the  analytical  result,  we 
can  solve  for  the  required  ratio.  Mote  aX  above  =  2aX$  jeon's. 
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This  ratio  is  applied  at  the  endpoint  of  the  integrals  in  Eqs.  (A-21)  and 
(A-22)  or  (A-29)  and  (A-30)  when  (Z-Z^)  >  fy-y^f,  there  is  a  lateral 
discontinuity  at  the  panel  edge  defining  the  lateral  limit  of  integration 
opposite  the  characteristic,  and  Y1  is  not  zero.  When  Y1  is  zero,  it  can  be 
shown  that  the  correct  limiting  value  will  be  obtained  by  using  half  of  the 
usual  non-integral  term  and  that  no  ratio  is  needed  on  the  integral  term. 

The  wing  solution  included  in  the  current  version  of  the  routine  involves 
some  approximations.  The  calculations  on  the  wing  utilize  a  routine  which  was 
developed  earlier  for  wing-alone  calculations.  It  uses  a  planar  source  sheet 
to  represent  the  wing,  and  the  boundary  conditions  are  satisfied  in  the  chord 
plane  rather  than  on  the  upper  and  lower  surfaces.  Therefore,  the  wing 
solution  is  symmetrical ;  i.e.,  it  is  the  same  on  the  upper  and  lower 
surfaces.  Only  the  normal  and  streamwise  perturbation  velocities  are  computed 
on  the  wing.  The  lateral  velocity  is  assumed  to  vary  from  zero  at  the 
wing-body  intersection  to  a  maximum.  The  maximum  occurs  at  the  leading  eoge 
for  subsonic  leading  edges  or  at  and  ahead  of  the  mach  line  from  the  wing-ooqy 
intersection  for  supersonic  leading  edges.  The  maximum  lateral  component  is 
such  that  the  total  perturbation  component  in  the  wing  plane  is  perpendicular 
to  the  leading  edge.  In  the  wing  solution,  except  for  the  effect  of  body  on 
wing,  a  reflection  plane  is  assumed  at  the  wing -body  intersection,  because  of 
the  planar  solution,  the  normal  velocity  is  equal  to  the  source  strength.  The 
wing -on-wing  solution  uses  the  same  modifications  to  linear  theory  as  the  body 
solution,  except  it  uses  see-through  rather  than  characteristic  tracing;  and 
the  free-stream  (X)  direction  integration  was  done  analytical ly .  For  the 
planar  case,  the  analytical  X  integration  case  simplifies  considerably  and  is 
more  efficient  computationally.  The  wing  slopes  and  source  strengths  are 
fitted  with  cubic  polynomials  streamwise  and  parabolas  spanwise. 

The  wing-body  Interferences  are  accounted  for  by  first  computing  the 
wing -alone  solution  with  a  reflection  plane  at  the  wing -body  intersection. 

Then  the  body  solution  is  calculated  including  the  effect  of  the  wing. 

Finally,  the  streamwise  perturbation  velocities  on  the  wing  due  to  the  body 
are  determined,  and  the  wing  solution  is  iterated  again  with  the  body-on-wing 
velocities  added. 

For  blunt  noses,  local  subsonic  flow  occurs.  The  modified  linear  theory 
will  provide  a  solution  with  local  subsonic  flow,  but  unless  the  local  Mach 
number  is  close  to  1.0,  the  accuracy  will  be  poor.  For  this  reason,  a 
modified  Newtonian  solution  is  included  in  the  method.  The  modified  Newtonian 
solution  Includes  an  adjustment  factor  from  Ref.  2  which  provides  improved 
accuracy  at  low  supersonic  freestream  Mach  numbers  on  hemispherical  noses. 
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The  pressure  coefficient  equation  Js: 

CB  = 


-  C.,  ( I  -  - '-<•« 

1  K'  3- r_  O.  3 


where  Cpn  is  the  pressure  coefficient  behind  a  normal  shock  wave.  M«,  is  the 
freestream  Mach  number,  and  s  is  the  local  slope  angle.  Hie  modified 
Newtonian  pressures  are  used  until  the  local  Mach  number  reaches  1.1.  Then,  a 
total  pressure  loss  factor  is  applied  to  match  the  modified  linear  theory 
pressure  to  the  modified  Newtonian  value.  This  pressure  loss  factor  is  then 
attenuated  exponentially  downstream.  The  scaling  constants  involved  in  the 
attenuation  were  determined  empirically  based  on  a  small  amount  of  data,  and 
further  analysis  based  on  a  broader  data  base  would  be  highly  desirable.  To 
prevent  the  (l-8e)  term  in  the  denominator  of  the  non- integral  terms  in 
and  ?2  from  reaching  zero  at  Be=l,  it  is  replaced  for  Be  >  0.8  with  an 
arbitrary  function  with  a  minimum  value  of  0.1  at  Be  =  1.0.  This  is  necessary 
to  provide  a  reasonable  solution  on  blunt  noses.  At  longitudinal  slope 
discontinuities,  the  exact  isentropic  solution  is  easily  calculated  in 
addition  to  the  modified  linear  theory  solution.  When  the  modified  linear 
theory  solution  predicts  a  higher  pressure  immediately  behind  the  corner  than 
the  two-dimensional  solution,  a  pressure  loss  factor  is  applied  in  the  same 
manner  as  was  done  at  the  modified  Newtonian  matching  point  above.  The 
pressure  loss  factor  is  again  attenuated  exponentially  downstream. 

The  values  of  source  strength,  pressure  coefficient,  velocity  components 
and  a  are  all  assumed  to  vary  linearly  across  a  panel  in  the  lateral 
direction.  They  are  continuous  across  panel  boundaries  unless  there  is  a 
lateral  slope  discontinuity.  In  any  given  row,  the  numDer  of  unknown  values 
of  source  strength  will  be  the  number  of  panels  plus  one,  plus  the  number  of 
lateral  slope  discontinuities,  tach  slope  discontinuity  adds  an  unknown 
because  at  that  panel  edge  there  are  two  values  of  source  strength  instead  of 
one.  The  number  of  control  points  at  which  the  boundary  condition  of  no-flow¬ 
through  is  satisfied  must  be  equal  to  the  number  of  unknowns.  Starting  at  the 
top  centerline  of  the  body,  control  points  are  placed  near  the  initial  panel 
edge  and  also  near  the  opposite  edge  of  only  the  last  panel  and  the  panels 
preceding  a  lateral  slope  discontinuity.  The  control  points  are  offset  from 
the  edge  by  a  fraction  (currently  equal  to  0.2)  of  the  panel  width  to  provide 
better  averaging  in  the  event  of  non-linear  variations  and  to  avoid  possible 
singularity  problems.  The  assumed  values  of  source  strength  and  b  are  at  the 
panel  edges.  The  calculated  values  of  perturbation  velocities,  B  and  pressure 
coefficient  are  at  the  control  points,  and  these  are  extrapolated  linearly  to 
the  panel  edges  for  interpolation  purposes  and  to  provide  the  8  estimates  for 
the  next  iteration. 


A  restart  capability  is  incorporated  to  allow  extra  prints  or  geometry 
changes  at  minimum  cost.  Using  the  restart  capability,  geometry  changes  can 
be  made  without  having  to  recalculate  the  solution  ahead  of  the  point  where 
the  geometry  changes.  After  the  calculations  for  each  row  of  panels  has  been 
completed,  the  pertinent  common  data  is  written  to  a  file  which  can  be  saved, 
providing  capability  to  restart  the  calculation  for  one  Mach  number  at  the 
beginning  of  any  row  for  which  all  preceding  row  computations  were  completed. 
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5*1  0  COMPARISONS  WITH  WIND  TUNNEL  DATA  AND  OTHER  THEORIES 


Five  cases  were  selected  for  testing  the  predictions  of  the  modified 
linear  theory  against  wind .tunnel  data  or  other  theoretical  methods.  The 
results  of  these  comparisons  could  then  be  used  to  establish  the  validity  of 
the  modified  linear  theory  method.  Comparisons  for  the  five  cases  are 
discussed  below. 

5.1  Analytic  Fore body 

An  analytically  defined  forebody  shape  was  tested  in  the  Langley  Unitary 
Plan  wind  tunnel.  Pressure  data  were  obtained  at  several  Mach  numbers, 
angles-of-attack,  angles  of  yaw,  and  at  many  body  locations.  These  pressure 
data  are  reported  in  Reference  3.  The  modified  linear  theory  method  was  used 
to  calculate  pressures  on  this  body,  called  Forebody  4  in  Reference  3,  at  zero 
angl  e-of-attack  and  zero  yaw  at  the  four  tested  Mach  numbers.  The  body  shape 
is  shown  in  Figure  7.  Calculated  pressures  are  compared  with  the  measured 
pressures  of  Reference  3  in  Figures  8  through  17.  The  agreement  is  quite  good 
considering  the  high  slopes  on  portions  of  the  body.  Linear  theory 
significantly  underpredicts  the  pressures  on  the  upper,  forward  portions  of 
the  forebody  at  Mach  1.7,  but  the  modified  linear  theory  compares  well.  At 
the  higher  Mach  numbers,  linear  theory  is  unusable  because  3e  is  greater  than 
1.0,  but  the  modified  theory  continues  to  compare  well  with  the  data. 

5.2  Drooped  Nose  Bodies 

A  series  of  bodies  with  circular  cross-sections  and  parametrically  varying 
amounts  of  nose  droop  have  been  tested  in  the  NASA/Langley  Research  Center 
Unitary  wind  tunnel.  The  results  are  to  be  published  in  the  near  future  by 
3arry  Shrout.  One  of  the  body  shapes  with  sixteen  degrees  droop  is  shown  in 
Figure  13  as  displayed  by  the  graphics  routine.  The  pressure,  or  wave,  drag 
results  were  obtained  by  subtracting  estimated  skin  friction  drag  from  the 
measured  forces.  Pressure  drags  predicted  by  the  modified  linear  theory  for 
three  bodies  at  three  Mach  numbers  are  compared  with  the  test  data  in  Figure 
19.  The  agreement  is  quite  good,  and  is  much  better  than  the  far  field  (area 
rule)  results. 

5.3  Conical  Wing  Body 

The  conical  wing  body  consists  of  a  symmetrical  wing  with  blunt  trailing 
edges  and  a  mostly  conical  body  underneath.  The  maximum  thickness  of  the  wing 
is  at  or  near  the  trailing  edge.  The  planfom  is  a  570  swept  delta  with  the 
outer  trailing  edge  cut  off  at  an  angle.  The  body  is  a  segment  of  an  8° 
cone  over  more  than  half  its  length  with  a  transition  approaching  a  constant 
section  near  the  aft  end.  Sec  Figure  20.  This  configuration  has  been  tested 
in  the  NASA/Langley  Research  Center  Unitary  wind  tunnel  and  the  data  are 
reported  in  Reference  4.  At  Mach  1.6  and  a  Reynolds  number  of  2.0  x  10 
*Yft. ,  the  Com-n  is  0.0186.  Subtracting  estimated  skin  friction  drag  of 
0.0063  leaves  a  measured  pressure  drag  of  0.0123.  The  modified  linear  theory 
calculation  yields  a  pressure  drag  of  0.0039,  consisting  of  0.0014  body  drag 
and  0.0025  wing  drag.  A  close  examination  of  the  geometry  reveals  that  the 
wing  leading  edge  is  very  blunt.  Also,  the  wing  is  conical  to  the  streamv/ise 
station  of  the  wing  tip.  A  conical  flow  calculation  using  the  modified  linear 
theory  with  the  wing  defined  as  part  of  the  body  and  input  points  concentrated 
near  the  wing  leading  edge  yields  a  drag  prediction  of  0.0082.  Thus,  most  of 
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the  discrepancy  is  explained  by  the  inability  of  the  wing  solution  to 
adequately  represent  a  shape  where  most  of  the  drag  occurs  in  the  first  one  or 
two  percent  of  the  chord. 

5.4  Missile  With  Inlets 


A  missile  model  with  axisymnetric  inlets  was  tested  at  NASA/Langley 
Research  Center  and  the  data  are  to  be  published  soon  by  Clyde  Hayes.  The 
arrangement  is  shown  in  Figure  21.  The  ,0=90°  tail  orientation  with  wings 
off  is  compared  at  Mach  2.5  with  the  wind  tunnel  data.  The  measured  Cq  . 
is  0.52  from  which  estimated  skin  friction  drag  of  0.20,  spillage  drag  8?n 
0.04,  and  boundary  layer  diverter  drag  of  0.02  are  subtracted,  for  a  net 
experimental  wave  drag  of  0.26.  Due  to  difficulties  with  tiie  computer 
routine,  it  was  necessary  to  run  the  configuration  in  parts:  the  forward 
body,  the  aft  body,  tails  and  inlet  lip.  The  modified  linear  theory 
prediction  is  0.27  for  the  wave  drag  on  the  body,  tails,  and  inlet  lip.  This 
is  considered  to  be  good  agreement. 

5.5  Area  Rule  Comparisons 

The  axi symmetric  body  for  this  case  consists  of  a  tangent  ogive  nose 
followed  by  an  arbitrarily  "coked"  section.  A  trapezoidal  wing  with 
leading  edge  sweep  and  a  65AQG5  airfoil  is  located  at  two  positions.  See 
Figure  22.  The  forward  position  provides  positive  interference  drag  between 
wing  and  body  and  the  aft  position  results  in  negative,  or  favorable, 
i nterference.  Area  rule  and  modified  linear  theory  drags  were  calculated  for 
Mach  numbers  from  1.2  to  2.0.  The  body  alone  and  wing  alone  drags  are 
compared  in  Figure  23.  The  body  alone  drags  agree  quite  closely  with  the 
modified  linear  theory  values  being  slightly  higher  at  Mach  1.3  and  2.0.  The 
wing  alone  drags  agree  fairly  well  except  at  the  sonic  edge  Mach  condition 
where  the  area  rule  drags  are  always  high.  At  Mach  1.0  and  2.0  (supersonic 
leading  edge),  the  modified  linear  theory  drags  are  slightly  higher. 
Interference  drags  and  total  drags  for  the  two  wing  positions  are  shown  in 
Figure  24.  The  interference  drag  predictions  are  similar,  except  the  modified 
linear  theory  does  not  predict  the  very  large  interference  effects  at  Mach  1.2 
as  does  the  area  rule.  Intuitively,  the  area  rule  interference  effects  at  1.2 
Mach,  particularly  for  the  aft  mounted  wing,  seem  too  large.  Also,  on  the 
total  drag,  the  modified  linear  theory  curves  are  smoother  because  of  the  more 
reasonable  results  for  the  wing  and  interference  drags. 

6.0  CONCLUSIONS 

Comparisons  between  v/ave  drag  and  pressures  as  predicted  by  the  modified 
linear  theory  and  as  obtained  from  wind  tunnel  tests  or  area  rule  predictions 
were  made.  In  the  conical  wing  body  case,  agreement  was  poor  because  of  the 
inability  of  the  wing  portion  of  the  routine  to  adequately  model  the  extremely 
blunt  leading  edge.  In  the  other  four  cases,  agreement  is  good.  For  the 
analytic  forebody,  good  agreement  was  obtained  with  wind  tunnel  measured 
pressures  over  a  wide  Mach  number  range  on  a  shape  which  has  high  local 
slopes,  where  linear  theory  is  unusable  and  area  rule  is  highly  inaccurate. 

On  the  drooped  nose  bodies,  pressure  drags  predicted  by  the  modified  linear 
theory  agree  well  with  the  wind  tunnel  data  for  the  three  droop  angles  and 
three  Mach  numbers.  For  the  missile  with  inlets,  wave  drag  calculated  by  the 
modified  linear  theory  agreed  with  wind  tunnel  data  for  the  one  case 
examined.  In  the  area  rule  comparisons,  the  modified  linear  theory  drags 
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agreed  well  with  the  area  rule  values  except  in  a  few  areas  where  there  is 
reason  to  believe  that  the  modified  linear  theory  results  are  more  realistic 
than  the  area  rule  values.  It  is  concluded  that  the  accuracy  of  the  modified 
linear  theory  has  been  verified. 

It  is  recommended  that  the  modified  linear  theory  be  used  to  predict 
supersonic  pressures  and  wave  drag  on  any  bodies  and  wing-bodies  which  include 
areas  of  moderate  to  high  local  slopes,  or  which  have  any  areas  including  wing 
leading  edges  which  have  supersonic  edges  (be  >_  1).  It  is  also  recommended 
that  the  computer  routine  be  improved  to  make  it  more  efficient,  general  and 
rel iable. 
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Upper  Surface  Pressures,  Mach  1,7 


Pressures,  Mach  2.5 


Fnrface  Pressures,  Mach  3.°5 


Figure  11  -  l;:per  Surface  Pressures,  Mach  4,5 


Figure  12  -  Station  Pressures,  Mach 


Figure  13  -  Station  Pressures,  Mach  1.7 
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Figure  14  -  Station  Pressures,  Mach  1.7,  X  =  10.0 
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Figure  17  -  Station  Pressures,  Mach  4.5,  X  =  10.0 
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Figure  19  -  Drooped  Nose  Body  Drags 
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Figure  23  -  Area  Rule  Comparison  -  Wing  Alone  and  Body  Alone  Drags 
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APPENDIX  A 


DERIVATION  OF  PERTURBATION  VELOCITY  EQUATIONS 


Curved  (Circular  Arc)  Panel  with  Source  Strength  Laterally  Constant 


The  velocity  potential  is 


,  _  f f'  5  (  <.}  sL  < ,  cLs _ 

rr.-  y  /(x-  x.y1  $.}*■-  ?-( *  -i,)4 


where  <<  f  -  /<*^,  *  z/  on  the  surface  of  the  panel.  No  subscript 

or  a  Q  subscript  (below)  indicates  the  receiving  point,  and  the  subscript  1 
indicates  the  sending  point.  S (x i )  is  the  source  strength  at  xi.  The 
curved  panel  is  assumed  to  be  a  segment  of  a  circle;  thus  with  ds  =  ri  d  u, 
yi  =  ri  sin  u  and  i\  =  ri  cos  u  we  have 
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where  y  and  z  must  be  defined  relative  to  the  center  of  the  circular  arc  panel 
segment.  On  any  panel,  the  curvature  and  (y,z)  coordinates  of  the  arc  center 
can  vary  with  x.  The  angle  u  is  positive  clockwise  and  r^  is  defined  as 
positive  for  convex  panels  and  negative  for  concave  panels.  We  will  integrate 
from  left  to  right,  so  rj  d  p  will  be  positive  for  both  convex  anu  concave 
panels.  Now  if  we  let  y  =  o,  the  m  =  o  direction  is  defined  by  the  vector 
from  the  center  to  (o,z)  for  convex  panels,  and  the  opposite  direction  for 
concave  panels  (i.e.,  from  the  receiving  point  to  the  center  of  the  panel 
arc).  Also  define  rq  as  the  vector  length  from  the  circular  arc  center  to 
the  receiving  point,  but  with  the  same  sign  as  rj.  Now 
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In  Eq.  (A-4),  if  |q  £  1,  the  inner  integral  is  an  elliptic  integral  of 
the  first  kind  with  modulus 


and  kj  <.  1  when  xi  £  x*,  where  x*  is  defined  as  the  xi  where  the 
forecone  passes  through  the  opposite  side  of  the  (imaginary)  complete  circle 
defined  by  the  radius  of  the  panel  element  being  considered,  or 


x*  =  x 
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For  x  >  x*,  we  must  transform  again  using  sin  u  =  kj  siny/  .  This  gives: 
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For  the  last  transformation  to  be  s  1  no le -valued,  we  must  limit  -  -t/2  <■„•<_  t/2 
(  -  ir  <  \i  <  n).  Then  cos  o  =  J  1-  sin^  u  and  cos  v  =  /  1  -  U..  iZ  ^  *- 


Defining  k£  =  1 /k i ,  we  have: 
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where  to  include  the  complete  area  of  integration,  the  bracketed  terms  in  Eq. 
(A-8)  must  be  summed  over  all  the  panels  within  the  area  of  integration 
defined  by  the  character) sties.  Also,  except  for  bodies  of  revolution,  x* 
will  be  different  for  each  panel. 

The  next  step  is  to  differentiate  with  respect  to  x,  y  and  z  to  obtain  the 
perturbation  velocities.  The  y  derivative  will  be  taken  first  because  it  is 


simpler  and  Is  the  same  for  xj  <  x*  and  xi  >  x*.  Going  back  to  Eq. 

(A-4),  the  only  terms  that  are  functions  of  y  are  the  limits  of  integration 
on  v  when  not  on  a  characteristic,  and  then 

y  2±. _ 


The  y  derivative  of  Eq.  (A-4)  can  be  written  in  several  forms: 
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but  if  either  lateral  limit  is  on  a  characteri stic ,  set  that  term  to  zero 
because  on  a  characteristic 

JL*'  r  P-1:'  =  o 

For  the  x  and  z  derivatives,  we  will  need  to  take  derivatives  of  the 
limits  of  integration.  The  upper  limit  (xr  or  xi)  of  x  integration  is,  in 
the  vicinity  of  a  point  A  on  the  surface  of  the  body  where  e  is  the  local  slope, 
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The  lateral  limits,  when  x  <  x*.  are  functions  of  y  or  which  are  independent 
of  x  and  z,  so 
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When  x  >  x*,  the  «  transformation  is  such  that  on  the  forecone 
(characteristic),  u  =  t  */2.  Thus, 
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But  with  x  >  x*  and  the  lateral  limit  not  on  a  characteristic,  from  sin  u  * 
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How  we  are  ready  to  differentiate  Eq.  (A-8)  in  three  parts:  11)  the 
derivative  of  the  upper  limit  on  the  x  integration,  (2)  the  derivative  of  the 
lateral  limits  of  integration,  and  (3)  the  derivative  of  the  integrand,  a 
contribution  to  part  (1)  occurs  at  xc  =  x  when  the  receiving  point  is  on  the 
surface  of  the  body,  and  also  whenever  there  is  a  local  maximum  of  the 
integration  variable  X}  due  to  the  forecone  from  the  receiving  point  tracing 
directly  to  the  body  (panel)  at  a  right  angle.  At  these  points,  =  u  and 
F(0,  1  ir/2)  -  F(0,  -  w/2)  =  z  n.with  the  sign  determined  by  the  sign  of 
ri.  So  differentiating  the  upper  limit  of  Eq.  (A-8)  and  substituting  Eqs. 
(A-12) , 
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where  the  forecone  from  the  receiving  point  traces  directly  to  without 
first  passing  through  or  around  a  body  and  the  area  inside  the  characteristic 
has  a  local  x  maximum  at  xj_.  The  z  direction  in  this  case  is  the  outward 
normal  from  the  sending  panel  at  xj.  In  the  most  cornnon  case,  xj  =  x  and 
rl  =  rQ. 

The  second  part  of  the  x  and  z  derivatives  comes  from  the  derivative  of 
the  lateral  limits  of  integration.  This  contribution  exists  only  when  x  >  x* 
and  the  limit  is  not  a  character!-  stic.  Differentiating  the  u  limits  of  Eq. 


Substituting  Eqs.  (A-15),  (A-16)  and  (A-5),  we  get 
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except  use  zero  when  q_  or  are  on  a  characteristic. 

The  third  part  of  the  x  and  z  derivatives  comes  from  differentiating  the 
integrand  of  Eq.  (A-8):  _ 
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Combining  the  contributions  from  Eqs.  (A-17),  (A-19  and  (m-20),  we  get  for 
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The  and  dz  integrands  are  functions  of  incomplete  elliptic  integrals  of 
the  first  kind  and  their  derivatives  with  respect  to  the  square  of  the 
modulus.  The  amplitude  may  be  negative,  and  F(k,  -  V  )  =  -  F(k.j').  The  dy 
equation  (A-10)  is  much  simpler  and  contains  no  elliptic  integrals. 

Several  pertinent  facts  regarding  Eqs.  (A-10),  (A-21)  and  (A-22)  are:  (1) 
the  non-integral  terms  occur  whenever  the  characteristic  from  the  receiving 
point  intersects  the  sending  panel  at  a  right  angle  without  around  or  passing 
through  any  intervening  surface  and  either  the  sending  panel  is  convex,  or  the 
sending  panel  is  concave  and  the  (y,z)  projection  of  the  receiving  point  lies 
between  the  panel  arc  and  the  center  of  the  arc;  (2)  the  integral  ten, is  are 
actually  lateral  summations  over  all  of  the  panels  within  the  bounds  of  the 
characteristics  and  X*  may  be  different  for  each  panel;  (3)  the  tenns 
containing  cos  u  must  be  set  to  zero  when  the  u  limit  is  a  characteristic;  (4) 
the  positive  z  and  u  =  ^  =  u  =  0  directions  are  defined  by  the  vector  from 
the  (y,  z)  center  of  the  sending  panel  arc  to  the  (y,  z)  projection  of  the 
receiving  point  if  the  sending  panel  is  convex,  and  the  reverse  if  the  panel 
is  concave;  (5)  the  moduli  of  the  elliptic  integrals  are  defined  by 


and  (6)  the  amplitudes  are  related  to  the  real  angle  u  by 


If'  =  1/2  u 

u  =  sin-1  ( k i  sin ) 
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Flat  panel,  as  used  here,  refers  only  to  the  lateral  shape.  The 
panel  may  have  longitudinal  curvature  and  twist.  The  positive  z  direction  is 
the  outward  normal  direction  from  the  sending  panel,  and  will  vary  with  x 
along  the  panel.  If  we  define  the  source  strength  as 


s^xl)  =  e0  +  ei  y 


(A  -  2J 


where  e0  and  ei  are  functions  of  x\,  the  velocity  potential  is: 


Carrying  out  the  lateral  integration  from  the  left  side,  y^_,  to  the  right 
side,  yR,  of  a  panel  we  get 


x . 


Now  we  will  differentiate  with  respect  to  x,  y  and  z  to  obtain  the 
perturbation  velocities.  Wherever  there  is  a  local  maximum  of  xc  or  xg, 
the  derivatives  of  the  limit  are  as  derived  in  Eq.  (A-12)  and  again  ~e>  ■< ,  : 

At  a  local  maximum  of  xc  or  x,  the  integrand  at  yR  minus  the  integrand  at 
yL  is  -  e0  it,  so  the  derivative  of  limit  terms  are 


When  yR  and  yL  are  not  characteristics,  they  are  not  functions  of  x,  y,  or 
z  so  there  is  no  derivative  of  lateral  limit  contribution..  When  yR  or 
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yi  are  on  a  characteristic,  the  integrand  is  infinite,  so  in  order  to  take 
only  the  finite  part  we  again  neglect  the  derivative  of  lateral  limit  term. 


Adding  the  x  derivative  of  the  integrand  to  the  derivative  of  limit  term 
in  Eq.  (A -2 8) , 


and  similarly  for  the  z  derivative, 


Note  that  the  integrand  is  identical  to  the  dx  integrand,  except  the 
sign  is  changed  and  (x-x^/gS  is  replaced  by  ( z-zi ) /s . 

The  only  contribution  to  tfy  is  the  derivative  of  the  integrand: 
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Several  pertinent  facts  regarding  Eqs.  (A-29),  (A-30)  and  (A-31)  are  : 

(1)  the  non-integral  terms  occur  whenever  the  characteristic  from  the 
receiving  point  intersects  the  sending  panel  at  a  right  angle  without  wrapping 
around  or  passing  through  any  intervening  surface;  (2)  the  integral  terms  are 
actually  lateral  summations  over  all  of  the  panels  within  the  bounds  of  the 
characteristics;  (3)  the  terms  containing 


v'  x  must  be  set  to  zero  when  the  y  limit  is  a 

characteristic;  and  (4)  the  positive  z  direction  is  defined  by  the  outward 
normal  from  the  sending  panel. 
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1.0  INTRODUCTION 


This  report  describes  Vought's  Improved  Wave  Drag  "Analysis  Routine",  ana 
an  interactive  computer  "graphics  routine"  to  view  ana  edit  the  input  aata 
deck  to  be  used  with  the  analysis  routine.  The  analysis  routine  computes 
local  velocities  ana  pressures  on  the  surface  of  the  configuration  ana 
integrates  the  pressures  to  obtain  lif-,  drag  and  pitching  moment.  At  the 
present  stage  of  development,  the  analysis  routine  can  hanale  a  single  ooay 
with  inlet  or  exhaust  surfaces  approximately  perpendicular  to  the  free  stream, 
and  one  pair  of  wing  panels  in  a  horizontal  plane.  The  winy  data  can  not  be 
input  to  or  displayed  by  the  graphics  routine.  The  data  deck  describes  the 
external  skin  contour  of  an  aircraft  fuselage  which  can  have  longitudinal 
steps,  and  may  have  a  blunt  nose  and  discontinuous  body  slopes  or  curvature. 
The  input  data  consists  of  three  dimensional  point  definitions,  together  witn 
sufficient  information  to  divide  the  points  into  rows  and  body  sections,  to 
identify  lines  of  slope  or  curvature  discontinuity,  ana  to  indicate  any 
"covered"  areas  where  no  boundary  conditions  are  applied  and  the  source 
strengtn  is  zero.  No  surface  slope  information  is  input  by  the  user. 

The  program  uses  a  tnree  dimensional  parametric  cuoic  space  spline  to 
define  slopes  along  curves  connecting  input  points  Doth  laterally  and 
longitudinally.  These  space  curves  then  are  used  to  define  "Coon's  patches" 
(parametric  cubic  surface  patches).  The  resulting  surface  definition  is 
guaranteed  to  be  continuous  in  position  as  well  as  first  ana  second  derivative 
except  along  lines  that  the  user  has  identified  as  slope  or  curvature 
discontinuities.  The  configuration  is  assumed  to  be  left-right  symmetric. 

Local  velocities  and  pressures  are  calculated  to  two  to  six  points  on  each 
panel,  and  the  source  strength  varies  linearly  in  both  directions  on  the 
panel.  The  source  strength  is  continuous  across  panel  boundaries,  unless  the 
surface  slope  is  discontinuous.  Because  of  these  techniques  and  the  curved 
panel  definitions,  the  panels  can  be  quite  large  in  regions  where  the 
velocities  and  source  stengths  are  not  expected  to  vary  rapidly.  For  example, 
in  a  region  where  the  configuration  is  nearly  ax i symmetric,  the  panels  coula 
span  thirty  degrees  or  more  laterally.  In  general,  the  input  deck  should 
contain  a  much  smaller  number  of  points  than  would  be  required  for  an  area 
rule  input  deck  describing  the  same  configuration.  A  permanent  restriction  on 


the  longitudinal  location  of  points  within  a  row  is  that  the  slope  (after 
rotation)  of  the  fore  and  aft  panel  boundaries  must  be  less  than  ,/Mj"l  / 

/I  +  t2  relative  to  the  Y-Z  plane,  where  Mf  iS  the  larger  of  true  local 
Mach  number  and  1.1,  and  e  is  the  panel  surface  slope.  The  program  will 
automatically  rotate  the  input  geometry  to  an  input  value  of  angle  of  attack. 

The  graphics  routine  has  been  implemented  on  a  CDC  6600,  CDC  Cyber  175, 
and  a  PRIME  mini-computer  with  a  Tektronix  4014  using  the  Plot  10  software 
package.  However,  the  program  was  written  to  be  easily  transportable  to  any 
type  of  graphics  scope. 


2.0  GEOMETRY  CONCEPT 

The  fuselage  external  surface  is  defined  by  a  set  of  three  dimensional 
point  definitions.  A  right-handed  X-Y-Z  coordinate  system  is  used:  positive  X 
is  aft,  positive  Y  is  the  right  side,  and  positive  Z  is  up.  Points  are  input 
in  rows,  and  the  resulting  set  is  connected  both  laterally  and 
longitudinally.  While  the  analysis  routine  only  requires  lateral  point 
connectivity,  the  splining  technique  used  to  determine  body  slopes  requires 
the  connectivity  in  two  directions.  Thus,  each  row  within  a  section  must  have 
the  same  number  of  points.  Points  are  ordered  from  bottom  to  top,  and  the 
rows  are  ordered  from  front  to  rear.  Lateral  symmetry  is  assumed  and  only  the 
positive  Y  half  of  the  fuselage  is  defined. 

The  wing  geometry  is  planar,  and  may  be  defined  by  airfoil  ordinates, 
standard  airfoil  definitions  or  streamwise  slope  inputs. 

3.0  INPUT  FORMAT  * 

3.1  INPUT  DATA  DECK,  EXCLUDING  WING  DATA 


The  data  deck,  except  for  the  wing  data  if  used,  consists  of  only  five 
types  of  cards:  title  cards,  point  definition  cards,  an  ENU  card,  an 
SrefCard,  and  a  Mach  number  card.  A  case  begins  with  two  title  cards.  The 
first  40  columns  of  both  cards  are  displayed  on  the  scope.  Following  the 
title  cards,  each  card  contains  a  single  point  definition,  with  five  free 
format  data  fields  in  the  first  40  columns.  Columns  41-80  can  be  used  for 
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comments.  Data  fields  are  separated  by  blanks  or  commas  ana  may  be  numeric  or 
alphabetic.  A  field  is  interpreted  as  numeric  if  the  first  character  is 
numeric  (0-9,  plus,  minus,  or  decimal)  and  is  alphabetic  if  the  first 
character  is  A-Z.  A  decimal  is  assumed  to  the  right  of  a  numeric  fiela  if 
none  is  present  within  the  field. 

The  five  data  fields  are  the  X,  V,  and  Z  values  for  the  point  and  two 
point  type  fields.  The  coordinate  fields  may  be  numeric  or  may  be  defined  by 
the  letters  "X“,  "Y"  and  “l",  to  indicate  to  use  the  previously  defined 
value.  The  point  type  serves  three  purposes.  It  indicates  the  division  of 
points  into  rows  and  sections,  indicates  slope  or  curvature  discontinuities, 
and  indicates  edges  of  covered  panels.  Point  type  is  an  alphabetic  field. 

The  first  character  indicates  whether  the  point  initiates  a  section  ("S"), 
initiates  a  row  ("R"),  is  a  corner  point  ("C"),  is  a  curvature  discontinuity 
("U"),  is  a  covered  area  ("V"),  or  is  a  normal  point  (anythiny  else;.  If  the 
first  character  is  an  "S",  "R"  or  "C",  the  second  and  third  characters  can 
indicate  the  direction  of  a  slope  discontinuity.  An  "X"  indicates  a 
longitudinal  discontinuity  and  a  "Y"  indicates  a  lateral  discontinuity.  Any 
other  characters  except  "D"  are  ignored.  For  example,  "CY"  indicates  a  point 
at  which  the  slope  is  discontinuous  in  the  lateral  direction,  and  "KX" 
indicates  the  beginning  of  a  row  at  which  the  longitudinal  slope  changes.  "C  " 

is  interpreted  as  “CY".  If  a  point  is  both  a  corner  point  (slope 

discontinuity)  and  a  curvature  discontinuity,  the  slope  discontinuity  must  be 
indicated  first  (e.g.,  CXDY  or  KYUX).  The  initial  point  of  a  data  set  is 
automatically  assigned  as  "S"  by  the  program.  If  less  than  four  fields  are 
supplied,  the  previous  values  are  assumed  for  the  missing  X,  Y  or  Z  fields  and 
"normal  point"  for  the  point  type.  Point  definition  input  is  terminated  by 
the  letters  "END"  in  the  first  data  field  of  a  card.  When  wing  data  are 

input,  they  are  inserted  after  the  body  data  are  completed  and  before  the  END 

card  (see  Sction  3.2  below).  Next  is  the  Sref  card  with  reference  area  and 
reference  chord  and  geometry  rotation  angle  (degrees)  in  three  data  fields 
(free  format).  There  must  be  at  least  two  non-zero  values  on  this  card.  The 
last  input  card  contains  analysis  Mach  numbers  in  free  format,  with  a  maximum 
of  19  values.  The  minimum  Mach  number  allowed  is  1.1  The  second  "Mach 
number"  may  be  used  as  an  indicator  -  e.g.,  U  for  only  wing-alone 
calculations,  -n  to  restart  at  row  n. 


The  covered  area  indicator  ("V")  appears  in  the  fifth  data  field  unless 
the  point  type  is  norami  (blank,  not  a  row  or  section  initiator).  If  the 
fourth  field  would  otherwise  be  blank,  the  V  goes  in  the  fourth  data  field. 

Any  additional  characters  in  the  V  field  are  ignored.  In  any  lateral  row  of 

points,  there  must  be  either  zero  or  two  V  indicators.  The  V  indicator  has  no 

effect  on  the  surface  splining.  When  the  V  indicators  are  matched  by  V's  on 
an  adjacent  row  at  the  longitudinally  connecting  pairs  of  points,  the  analysis 
routine  is  triggered  to  bypass  satisfying  the  no-flow-through  boundary 
condition  on  the  intervening  panel  or  panels  and  to  set  the  source  strength  on 
those  panels  to  zero. 

When  the  fuselage  cross-section  is  a  circle  and  seven  evenly  spaced 
lateral  points  will  match  up  with  adjacent  rows  in  the  section,  a  simple 

radius  and  Z  location  can  be  input  instead  of  the  seven  X,  V,  2  values,  (.I>«.» 

one  card  can  define  a  complete  row.)  In  this  case,  the  five  data  fields  are 
X,  radius,  "AX",  point  type  and  Z  center.  Tne  "AX"  in  the  third  field  may 
have  additional  characters  after  the  AX,  which  will  be  ignored.  The  point 
type  can  only  be  S,  SUX,  R,  RX  or  RDX.  The  Z  center  shifts  the  center  of  the 
section  vertically,  if  desired.  The  program  computes  the  X,  Y  and  L  of  the 
seven  evenly  spaced  (30  degrees)  lateral  points. 

3.2  WING  INPUT  DATA  DECK 


The  additional  input  data  for  the  wing  is  inserted  after  the  body  data  is 
completed  and  before  the  END  card.  The  wing  data  begins  with  WING  in  the 
first  free  format  data  field  of  a  card.  At  present,  the  wing  data  cannot  be 
included  in  the  input  to  the  graphics  routine. 

« 

When  wing  data  are  input,  the  rotation  angle  of  the  Sref  card  must  be 
zero,  and  there  must  be  only  one  real  (i  1.1)  Mach  number.  If  only  a 
wing-alone  calculation  is  desired,  input  two  Mach  numbers  with  the  second  one 
equal  to  zero.  The  reference  area  input  to  or  calculated  by  the  wing  routine 
will  override  the  reference  area  on  the  5pef  card  described  bove.  Also  note 
that  the  drag  values  printed  with  the  body  calculations  do  not  include  tne 
wing  drag  and  vice  versa. 

The  wing  surface  is  defined  by  airfoil  sections.  The  input  can  consist  of 
an  area  rule  deck  or  surface  slopes  which  will  be  curve-fitted.  Tne  wing 


planform  may  be  trapezoidal  or  it  may  have  one  leading  edge  break  (two  sweep 
angles)  and  up  to  three  trailing  edge  breaks  (four  trailing  edge  sweep 
angles).  The  wing  geometry  is  flattened  into  a  Z=constant  plane  through  the  Z 
value  of  the  root  chord.  There  is  an  upper  limit  of  six  spanwise  curve  fit 
intervals,  which  means  that  with  an  area  rule  input,  a  maximum  of  7  span 
stations  can  be  input  unless  a  discontinuity  flag  or  duplicate  span  station 
locations  (Y,  Z)  are  used. 

The  wing  inputs  are  described  in  the  following  figures  and  definitions. 
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WING  DATA  INPUTS 
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IF  FNORD  <3..  THE  X/C  CARD  IS 
OMIT-EC.  THE  PROGRAM  SETS 
F  N0R3  =  13-  AND  USES  EVERY 
10  •  CHORD. PLUS  5  7C  AND  15*£ 


1-3 

103 

6 

OR  7 


11 - 20 

NAME 

(ADDED  TO  END  OF  TITLE! 


A3, 


NOTE  Xle.Y.Z.  AND  CHORD  VALUES  MUST  BE  SUCH  THAT  THE 

RESULTING  [FLATTENED!  WING  HAS  NO  MORE  THAN  2  L-E. 
AND  4  T.E.  SWEEP  ANGLES. 

TOLERANCE  IS  A  (  T  AN.A  £  )c.0174S.A(TANArtX. 01745 


AREA  RULE  INPUTIS4007) 
DATA  ON  THIS  PAGE  IS  USED  ONLY  IF  IARULE=l 
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NOTE  NS  ARE  NEEOEO  ONLY  IF  AIRFOIL  HAS  CORNERS  AND  NWAF0R>5- 
IF  N’S  ARE  USED  THE  CHOROWISE  C'JRVEFIT  OF  THE  AIRFOIL 
SLOPES  HILL  BE  DISCONTINUOUS  AT  XAF.  ;  I  - E .  .  THE  VALUES 
OF  XAF  DEFINE  THE  CHCRDWISE  CURVEFIT  BOUNDARIES. 


NOTE  PROGRAH  SHOVES  LEFT  AND  RIGHT  WING  PANELS  TOGETHER 

if  this  is  nct  desired  but  rm>o.  2  additional  airfoils 

AT  YrO  AND  Y  =  Yf 1  )  WITH  ORD I NATESrO  HUS T  BE  ADDED. 

ALSO  PROGRAM  FLATTENS  The  WtNGI-iYr/aY.;^  *J.Z-  .  ) 

A - 7  S  — 14  15 — 21  22—23  C9--3S  35  —  42  43  —  45  SC— 56  S7--63  64  —  70  " 

■  WAFCRC,  Wfir  OhD.  WAFCRD,  WAFORD.,  WAFORD.  WAFORD.  WArCP.D  WAFORD  WAFORD  WAFORD. 
SEMI -THICKNESSES  Or  THIS  AIRFOIL 

CHIT  IF  STANDARD  AIRFOIL  OPTION  (QFID)  IS  USED 
WAFORD"1  S  CONTINUE  ON  SUBSEQUENT  CARDS  IF  NWAFOR  >10 


8 14  IS--2:  22  —  28  41 - SO  51 - 6C 

Y  Z  CHORD  AFIC  XTMAX 

OR 
T/C 

REPEAT  THIS  CARO  NWAF  TIMES 


61 - 70  73  7S  79 

CISC.FLG  N  N  N 


-7  8—14  IS— 21  22--2e  29- -35  36—42  43  —  49  50—56  57—63  64--70 


1QF7-0 

REPEAT 

NWAF 

TIMES 

1  4F7.0. 
i  12X. 
3F10.0. 
313 

|  10F7.D 


OMIT  IF  NWAFOR  <  3 

XAF’S  CONTINUE  ON  SUBSEQUENT  CARDS  IF  NWAFOR  >10 

fg  jo 

NWAF  NWAFOR 
(2  TO  12)  13  TO  25) 

i-BUT  FOR  NWAF  >7  TO  BE  ALLOWABLE. 

EITHER  DISC-  FLAG  OR  TWO  AIRFOIL  QEFS.  AT  THE 
SAME  V.Z  MUST  BE  USED  SO  THAT  6  OR  LESS 
SPANWISE  CURVEFIT  INTERVALS  RESULT. 

NOTE  X„  .Y.Z.  ANO  CHORD  VALUES  MUST  BE  SUCH  THAT  THE 

RESULTING  I  FLATTENED )  WING  HAS  NO  MORE  THAN  2  L-E. 
AND  4  T.E.  SWEEP  ANGLES. 

TOLERANCE  IS  acTAN^-H*.  01745,  At  TAN-A-^X  .  01745 


INTEGER 


AREA  RULE  INPUT! TMX-947) 
DATA  ON  THIS  PAGE  IS  USED  ONLY  IF  IARULE*2 
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IF  JCONIC  >  0 ,  ix  IS  INTERPOLATED  AS  A  FUNCTION  OF 
X- ( E _  *F _  «Y) 

TET'F_  -V  )^-t e _ 7  *f _ _  -n 

INSTEAO  OF  X/C 

-10  1  1 . 20  21 - 30  31 . 40  41 - SO  51 - 


6F10-0 


3  3  s  s  OMIT  IF  NSPPLP  =  1  OR  2  \\ 

THIS  CARD  USED  ONLY  IF  NPN>2( ADD  I T I ONAL  CARDS  NOT  SHOWN  ARE  USED  IF  N(’N>5  I  \\ 

'1 - 10  1  1 - 20  21 - 30  31 - 40  A  1 - SO  51 - 60  |  6F10-0  \\ 

Y3  ALP-MULT  (X/C)  IX,  (X/C),  oc.  \\ 

'  '  -  1  OMIT  IF  NSP ^LP- 1  OR  2  \ \ 

NOTE;  DO  NOT  PUT  AN  INTEGER  IN  COLUMN  ltCOL-l  IS  PART  OF  ThE  Y3  FIELO) 

/I - 10  1  1 - 20  21 - 30  31 - 40  41 - 50  51 - 60  H  6F  1 0 . 0  \ 

3  (X/C),  CX  (X/C),  «.  (X/C)  ex.  \ 

J  *  J  J  onIT  j  r  MSPfiLPlt  \ 

THIS  CARO  USED  ONLY  IF  NPN>2( ADD  I T I ONAL  CARDS  NOT  SHOWN  ARE  USED  IF  N »N>S )  \ 

/; - 10  n - 20  21 - 30  31 - 40  41 - SO  SI - SO  ’  |  6F10.0 

/  Y2  ALP-MULT  (X/C)  o/  ( X/C )„  <y, 

f  ’  *  ‘  OMIT  I -  NSPALP=1 

!  NOTE;  00  NOT  PUT  AN  INTEGER  IN  COLUMN  1 ( COL  .  1  IS  PART  OF  THE  Y2  FIELD) 


'  1 - 10  11- 

!  X/C  ), 


-20  21 - 30  31 - 40  41- 


SF10.0 


this  carc  used  only  if  npn>2( additional  cards  not  shown  ape  used  if  n|=n>5 )  1 

A  2 - 10  It - 20  21 - 30  31 - 40  41  - 50  E  1  -----SO  '61  - - 70  "l  1  .  F  9 . 0  . 

.  Y 1  ALP-MULT  !  X/C )  «  (X/C)  ex  JCCNIC.  5F10.0.I10 

NPN  - - 1 - ,  L  1  ALTERNATE  i 

(1  TO  9)  I  NOTE  o<  IS  IN  RADIANS  I  CUE  .'EF  IJ  OPUCN _ 

1 - - -  xXrEAT  fJFlNT  TIMES 


3  4 - 1C  II - 20  21 - 30  31 - 4C  41 - 50  51 - SC  61 - 70  ll.I2.7X. 

.  9LANK  El  FI  E2  F2  E3  F3  6F10-0 

NSPALPlNCRMALLY  1  TO  3. IF  LEFT  BLANK  IS  SET  TO  3.  BU'  IF  A  NEGATIVE 
va:.'-E  ISVnput.  Then  t.  Curve-fit  coefficients  will  BE  READ  IN  6Fio.) 


F  ORM  A  r 
I  NT 
TO  4) 


INSTEAD  OF  X/C. «  CATA  ABOVE) 

N — IF  9LANX  IS  NOT  BLANK. THE  CURVEFIT  BOUNDARIES  ARE 

DEFINED  AS  X/C'S  , _ _ 

I  REPEAT  N-EPAN-INTS  TIMES 


6F10-0 


AY( I  1  IS  THE  LEFT-HAND  OR  INBOARD  BOUNDARY 
OF  THE  I  TH  CURVE-FIT  INTERVAL:  AN  «  IN  THE 
FIRST  COLUMN  OF  AM  AY  FIELD  (11. 21. ETC-) 
INDICATES  A  SPANWISE  SLOPE  D I SCCNT I NU I TY- 
RESULTS  IN  AN  EXTRA  Y  AND  SPECIAL  INTEGRATION 
TREATMENT 

THIS  CARO  OMITTED  IF  N-SPAN-INTS=1 


X.Y,Z  COORDINATES  OF  ROOT  CHORD  LEAOINO  EDGE 


ALL  DATA  ON  THIS  PAGE  OMITTED  IF  IARULEPO 


DEFINITIONS 


Notation 

Synbol  Definition 

N-SPAN-IMTS  Number  of  ep.mwise  curve-fit  intervals 

ITHICK  Thickness  option;  ITHIC-l  =  0  for  camber, 

twist  and  tuple  of  attack  effects, 

ITHICK  =  1  for  thickness  effects 

NEXTRY  Number  of  extra  span  stations  (in  addition 

to  those  determined  by  multiples  of  DEL-Y) 
for  aero  calculations 


No  KEF 


I.e'uranc e  area  option;  NC2.EF  =»  0  total  plan 
«’ i  -a  will  :••»  uaed;  N3RX"  ••*  1  input  value 
o2  SEE F  './ill  used  in  scare  units  as  plan- 
fora  direr,  ices;  N2XZ"  =  2  input  value  of 
E.'d'.F  will  re  used  assumir-;  SEE?  is  in  ft2 
d  dimer.-: :  - r.  a  are  in  inch.  -s;  NS  REF  «  3  input 
v.> lue  of  O used  as  a  ratio,  i.e. ,  reference 
cir-a  »  S'"  plan  are  .;  NSREF  3  4 

reinrcr.o.:'  a  roe  “  plan  area  plus  an  extension 
in'. card  of  SEE?  inches  (or  whatever  units 
are  be lap  rv-J) ,  this  is  useful  for  thickness 
drag  cases  wee  re  exposed  panels  have  been 
shoved  topntr  -r  for  reflection  plane  effect; 
NSREF  =  5  reference  area  defined  by  outer¬ 
most  win;;  panel  extended  in  to  the  centerline; 
MS REF  “  6  reference  are3  defined  by  innermost 
wins  panel  extended  out  to  the  tip  station. 


MODE IK  Mon  linearity  option;  pressure  coefficients 

calculated  from  the  exact  equation  usinp 
calculated  strearvise  perrubntion  velocitv 
and  ass u:no d  srar.-.-Lse  pertubation  velccitv. 

In  addition,  when  ITUI0E.=  l,  local  ft  (fore- 
core  share)  is  modified  to  account  for  local 
perturbed  velocity  and  the  flow  tan coney  con¬ 
dition  us>s  local  velocitv  rather  than  free- 
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red 
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-n 
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■pertni!  t*x  • 


Units 

None 

None 

None 


None 


None 


L'o 1.  tuition 


Units 


t 


.  .  M.U 

Ni-ESX  P3 


No  •!'•  ••  of  SVC  op  angles  used  to  define 
loo  din;*,  odja  of  pi  an  torn 


Number  of  swoop  angles  used  to  define 
trailing  edge  of  planform 


I.V’.ULX 


C'orocrv  and  slrpo  input  ontion.  If  I A  PULE 
>  0  (allowed  only  if  ITHICK  >0,  JSYM  =  0, 

^  0 ,  and  EONIC  =  0),  the  wing  defini¬ 
tion  portion  of  an  area  rule  innut  deck 
provides  both  the  wing  geore  try  and  airfoil 
(slope)  inputs.  If  IARULE  =  1,  Vougnt  format 
is  used  and  if  IARULE  =  2,  NASA  (HI  X-947) 
format  is  used.  If  IARULE  <  0,  slope  inputs 
are  required.  If  IARULE  <  0,  geometry  is 
defined  by  (Y,  X)  coordinates.  If  IAXULE  = 

0,  geometry  is  defined  by  lending  and  trailing 
edge  s*.  eep  angles 


YEX(I) 


Extra  span  stations  for  aero  calculations 


None 


None 


None 


i 
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Unt  ts 


Symbol 

XXL ( I ) 
YLE(I) 

XXT(J) 

YTE(J) 

SWL  (I) 

SWT(J) 

DEL-" 

y-al:-?isc 


Dif  i  ni  tiers 


X  coc-rdiante  of  rhc  leading  edge  at  YLE  (I) 

Spaa  -station  defining  leading  edge 
bread  point. 


X  coordinate  of  the  trailing  edge  at  YTE(J) 

Span  station  defining  trailing  edge  break 
point 

Leading  edge  sweep-back  angle, 

S'.vL(n)  is  the  sweep  inboard  of  YLE(n)  . 


Sweep  defining  trailing  edge  of  planforcs 
SWT(n)  is  sweep  inboard  of  YTE(n). 


Increment  in  y  between 
station.;  if  left  blank. 
0.1  -  :r:n,  i:  -.-ill 

0.1  ti-  r.  -----  ’  -  -- 


calculated  span 
or  made  less  than 
be  set  equal  to 


Sp.u.  rtati'n  at  the  constant  (fuselage) 

or  tr.f-r-.--s  1  ends;  - 1  be  oxnc“ly  ecuai  to 
A:(-,  i-. 


•k 

k 

k 

k 


Degrees 


Degrees 


k 


k 


'.'one 


DEL -f LORE 

Ccn.i:  ir\ z ,  co  s--  added  no  input: 

Radians 

see" 

See  ' 

- 

CRCCT 

Root  Chord. 

k 
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Ue li uir  s on 

l!ni  ts 

AY 

Sp.nwt  ;;o.  curve  fit  boundaries;  value  is 
equal  to  the  v  value  of  the  left  boundary 
of  the  respective  span  interval. 

* 

C 

Twist  case  only  (JTWIST  >  0  and 

K-SPAA- 1NTS  »  1),  cv  =  CyJ'f»ISI. 

Radians/ 
(*)  JTWIST 

nxint 

Number  of  chordvise  curve  fit  intervals 
in  a  given  spanvise  interval 

None 

ns palp 

Number  of  span  stations  for  cc  inputs 

None 

El,  FI 

Coefficients  of  boundary  between  the 
first  and  second  cnordwise  intervals 
where  X  =  El  +  FI  (y) 

* 

e:.,  f: 

Coefficients  of  boundary  between  the  second 
and  third  chornwi.se  intervals  where 

X  -  Ed  +  F2  (y) 

* 

r  “>  ir  n 

-  - 

Coefficients  of  boundary  between  the 
third  and  fourth  chordvise  intervals 
where  X  -  E3  +  F3  (y) 

Vf 

NTS 

Number  of  input  pairs  of  X/C,  cd  per 
span  scarier.  for  the  particular  curve- 
fit  Interval  bring  considered 

Hone 

VI  '  ">  v  "i 

lx,  1  . ,  : - 

Tnc  three  span  stations  at  which  ■:•£  data 
is  given  far  the  particular  curve-fit 
interval  being  ccr.uidarcd 

* 

ALT *-)!..  .. 

Factor  by  which  input  ex's  for  tins  Y  and 
this  chordvise  interval  only  will  be  multiplied, 
hot  used  if  input  value  =  0 

None 

x/c 

Fraction  of  local  chord 

None 

D<_ 

Local  angle  of  attack  (for  ITHICK  =  1, 

«•  =  -d(t/2)/d::>. 

Radians 

JCONIC 

Curvefit  option;  if  70,  slopes  will  be  fitted 
and  interpolated  snar.wise  along  constant  fractions 
of  curve-fit  panel  chord  instead  of  wing  chord. 

This  option  is  chosen  separately  for  each  curve- 
fit  interval. 

Whe;;  an  opt-cn  mricator  it;  0,  the  cutler.  is  not  used,  when  it  is  1,  the  option 
is  used . 

♦Indicates}  that  th-  length  dl .-■•r.aiar.r.  -ast  merely  be  consistent. 
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4.0  INTERACTIVE  PROGRAM  CONTROL  (GRAPHICS  ROUTINE) 


4.1  command  summary 


For  reasons  of  initial  program  portability,  all  program  interaction  is 
handled  through  the  keyboard.  ATI  commands  are  initiated  by  typing  at  least 
tne  rirst  two  characters  of  the  command  name.  The  routine  will  respond  with  a 
prompter  describing  the  additional  information  (if  any)  required  to  complete 
the  command.  Commands  fall  into  four  categories: 

•  Ui splay  options 

•  Editing  options 

•  Editing  aids 

•  Output  options 

The  display  options  give  the  user  control  over  the  appearance  of  the 
.••3Prm:al  n  splay.  Options  are  availaole  to  view  the  model  rrom  any  angle,  to 
■  the  scale  and  view  center  of  the  display,  to  clip  the  display  to  a 

v  r-.-e  : -mensiona  1  view  box,  to  display  either  or  Doth  sides  of  tne  model,  anu 
;■  tne  screen  and  re -display  the  model. 

iti  editing  options  allow  tne  user  to  change  tne  title,  move  one  or  more 
mts,  c hinge  tne  corner  cooes,  or  to  reject  the  effect  ot  the  previous 


editing  aids  allow  the  user  to  mark  a  number  of  points  on  the  display  with 
*h-.-i>-  y;iot  numbers,  or  to  display  the  coordinates  and  corner  codes  for  one  or 
mi. re  points. 

output  options  include  the  ability  to  print  or  punch  the  current  data  set, 
to  plot  the  display.  The  "plot"  option  currently  just  writes  a  vector  file 
wiicn  can  pe  processed  by  an  external  routine  to  route  the  information  to  an 
ic  p '-op  ri  ate  plotting  device.  Two  examples  of  the  scope  display  are  shown  in 

'idu'rs  1  and  d. 
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4.2 


DISPLAY  COMMANDS 


For  display  purposes,  the  complete  model  definition  is  stored  as  a  single 
vector  list,  in  model  coordinates.  When  the  program  calls  for  the  oi splay  of 
the  model,  each  vector  is  processed  by  the  following  transformations: 

1)  Translate  with  respect  to  a  view  center 

2)  Rotate  through  yaw,  pitch  and  roll  angles 

3)  Scale  to  a  specified  unit  sphere 

4)  Clip  to  a  specified  three  dimensional  viewport 

The  remaining  visible  portion  of  each  vector  (if  any)  is  displayed  with 
the  X  coordinate  mapping  onto  the  horizontal  screen  dimension  anu  the  L 
coordinate  mapping  onto  the  screen  vertical  dimension. 

The  display  options  allow  the  user  to  specify  the  parameters  used  in  the 
viewing  transformations  as  well  as  to  control  re-oi splaying  the  model  ano  to 
select  one  or  both  sides  of  the  model  for  display. 

Since  the  user  may  wish  to  change  several  of  the  display  parameters  before 
re-displaying  the  picture,  an  optional  display  suppression  flay  may  oe  input 
following  the  last  mentioned  parameter.  When  any  non-blank  value  is  supplied 
for  this  flag,  the  picture  will  not  be  re-displayea  until  a  further  command 
specifically  requires  it. 

4.2.1  Translate 

The  program  responds  with  "ENTER  VIEW  CENTER  X,  Y,  Z".  The  user  shoulo 
type  the  model  coordinates  of  the  point  that  will  be  translated  to  the  center 
of  the  screen. 

4.2.2  Rotate 

The  program  responds  with  "ENTER  YAW,  PITCH,  ROLL".  The  user  should  type 
the  values  in  degrees.  The  (0,  U,  0)  orientation  is  the  left  profile. 

Positive  yaw  is  nose  left.  Positive  pitch  is  nose  up.  Positive  roll  is  right 
side  up. 


A  convenient  way  to  visualize  a  desired  rotation  is  to  image  that  the 
viewer  is  initially  located  to  the  left  side.  First  move  around  towara  the 
nose  by  the  desired  yaw  angle.  Then  move  the  imaginary  viewpoint  through  a 
desired  elevation  angle.  The  pitch  angle  will  be  the  negative  of  the 
elevation  angle  and  roll  will  be  zero.  The  model  can  be  viewed  from  any  angle 
with  only  two  angles  specified.  Three  rotation  angles  are  necessary  only  to 
control  the  orientation  of  the  resulting  view. 

4.2.3  Scale 

The  program  responds  with  "ENTER  VIEWING  RADIUS".  The  user  enters  a 
radius,  in  model  coordinates,  that  will  be  scaled  to  fit  within  a  12-incn 
viewing  area. 

4.2.4  Clip 

The  program  responds  with  "ENTER  XMIN,  YMIN,  21-llN,  Xi-iAX,  YI-iAX,  ZI-iaX".  Tne 
values  are  in  display  coordinates  (screen  inches).  (U,  U,  0)  is  the  screen 
center.  The  following  may  aid  in  specifying  tne  values. 

XMIN  -  Left  side 

YMIN  -  "Near"  side.  Negative  Y  is  toward  the  viewer 
ZMIN  -  Bottom 
XMAa  -  Right  side 

YMAX  -  "Back"  side.  Positive  Y  is  behind  tne  screen 

ZMAX  -  top 

4.2.5  Side 

This  controls  whether  the  program  displays  the  left  side  l-Yj,  the  right 
side  (-*-Y ) ,  or  both  sides  of  the  model.  Tne  program  responds  "ItFT,  KlbHT  UR 
BOTH". 

4.2.6  Paint 

This  causes  the  program  to  clear  the  screen  ano  re-aisplay  the  picture. 
This  may  be  necessary  because  of  an  accumulation  of  command  or  output  data  on 
a  Tektronix  scope. 
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4.3 


EDITING  AIDS 


The  commands  are  available  to  aid  the  user  in  identifying  and  locating 
points  on  the  display. 

4.3.1  Hark 

This  command  causes  the  program  to  show  the  point  number  of  one  or  more 

points  on  the  display  at  the  point  location.  This  enables  the  user  to 
identify  points  in  the  configuration.  The  program  responds  with  "ENTER  POINT 
RANGE".  The  user  types  in  the  first  and  last  numbers  to  be  marked.  Once  a 
range  of  points  has  been  marked,  the  command  stays  in  effect  through  changes 
of  display  until  changed.  To  turn  off  the  point  marking,  enter  null  print 
range  (0,  0  or  only  a  space  before  the  end-of-1  ine)  . 

4.3.2  Locate 

The  program  responds  "ENTER  POINT  RANGE".  After  the  user  types  in  the 
first  and  last  point  of  interest,  the  program  will  display  the  existing  X,  Y 
and  Z  coordinates  and  the  corner  code  for  the  specified  points. 

4.4  EDITING  COMMANDS 


Four  commands  are  available  that  actually  change  the  data  set.  These 
allow  the  user  to  move  points,  change  the  corner  codes,  reject  a  previous 
change,  or  change  the  title. 

4.4.1  Hove  Points 

The  program  responds  "ENTER  POINT  RANGE,  X,  Y  Z".  The  user  enters  the 
first  and  last  point  to  be  moved  (which  may  coincide)  and  the  new  point 
coordinates  for  the  first  point  of  the  range.  All  the  points  will  then  be 
moved  by  the  same  increment  as  the  first  point.  This  allows  an  entire  row  or 
group  of  rows  to  be  moved  as  easily  as  a  single  point. 
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4.4.2  Change  Corners 


The  program  responds  "ENTER  POINT  RANGE  AND  NEW  CORNER  CODE".  The  user 
enters  the  first  and  last  point  of  interest  and  the  new  ten  character  point 
type,  which  can  include  the  regular  point  type  field  and/or  the  "V"  indicator 
for  covered  panels.  Notice  that  the  routine  will  not  redefine  rows  or 
sections  with  this  command,  only  the  slope  discontinuity  indicators  are 
changed. 

4.4.3  Reject 

This  command  retracts  the  effects  of  the  previous  HOVE  or  CHANGE  command. 

No  additional  parameters  are  entered. 

4.4.4  Title 

The  routine  responds  "TYPE  IN  FIRST  LINE  OF  NEW  TITLE"  and  then  "TYPE  IN 
SECOND  LINE  OF  NEW  TITLE".  Forty  characters  may  be  entered  in  each  line. 

4.5  OUTPUT  COMMANDS 


Three  output  options  are  available  to  allow  the  user  to  print  or  punch  the 
current  data  set,  or  to  write  a  plot  file. 

4.5.1  Print 

No  other  parameters  are  entered.  The  printed  output  represents  the 
current  title,  point  and  corner  codes. 

4.5.2  Punch 

No  other  parameters  are  entered.  The  punch  deck  may  be  used  for  direct 
input  either  to  the  graphics/editing  routine  or  to  the  analysis  routine  later. 
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4.5.3  Plot 


No  other  parameters  are  entered.  The  routine  writes  a  binary  file 
consisting  of  the  current  vector  set  and  view  transformation  parameters  for 
use  by  an  external  batch  routine  to  direct  the  required  plot  data  to  the 
appropriate  plotting  equipment.  If  the  vector  set  has  not  changed  since  the 
previous  plot  was  requested,  only  the  new  view  parameters  are  written. 

5.0  SAMPLE  CASE 

Printout  for  a  sample  case  is  shown  in  Figure  3.  Included  in  the  printout 
is  a  side  by  side  literal  and  interpreted  print  back  of  the  fuselage  input 
data.  If  a  wing  is  input,  the  print  back  is  the  interpreted  type  only.  Next 
is  the  panel  geometry  data  resulting  from  the  three  dimensional  splining. 
Finally,  the  computed  local  pressures  and  velocities  and  the  integrated  forces 
are  printed. 


6.0  SUbliEST IUliS  FUK  INPUT  PKEPAKATIUH 

The  most  important  recommendation  is  to  always  use  the  graphics  routine  to 
check  and  edit  the  body  geometry.  This  step  nearly  always  uncovers  one  or 
more  errors  in  the  input  data  or  areas  where  the  curve  fit  is  not  quite  what 
was  desired. 

In  general,  one  should  use  the  minimum  number  of  points  which  adequately 
describe  the  geometry.  The  only  exception  is  in  areas  where  it  is  apparent 
that  flow  conditions  will  be  varying  rapidly  due  to  asymmetry  or  local  slope 
variations.  In  areas  where  flow  conditions  are  nearly  constant,  a  panel  can 
span  30,  45,  or  even  90  degrees  laterally.  The  panels  should  not  span  more 
than  90  degrees  because  that  could  cause  inaccuracy  in  the  characteristic 
tracing. 

The  body  points  are  input  from  bottom  to  top,  starting  with  the  front 
now.  Although  the  input  points  are  sequenced  from  bottom  to  top,  the  panels 
are  numbered  from  top  to  bottom  and  the  results  are  printed  in  top  to  bottom 
order.  Each  lateral  spline  fit  extends  from  the  bottom  point  to  the  first 
lateral  corner  (CY)  or  the  top  point,  and  the  next  one  to  the  next  CY  or  the 
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top  point,  whichever  occurs  first.  Each  longitudinal  line  is  splined  from  the 
first  row  in  the  section  to  the  next  longitudinal  corner  (CX  or  RX)  or  to  the 
last  row  in  the  section,  etc.  A  curvature  discontinuity  (DX  or  DY)  does  not 
interrupt  the  fit;  it  keeps  the  slope  continuous  but  substitutes  a  constant 
curvature  condition  in  an  adjacent  segment  instead  of  the  equal  curvature 
condition  across  the  point.  Within  each  spline  fit  of  more  than  two  points 
there  must  be  at  least  one  point  at  which  the  curvature  is  continuous. 

The  discontinuous  curvature  (DX  or  DY)  option  should  be  used  when  one 
class  of  curve  joins  another  with  the  slope  continuous.  An  example  is  a 
tangent  ogive  cylinder  body  shape.  Curvature  reversals  also  can  occur  at 
locations  where  it  is  necessary  to  keep  the  slope  continuous. 

The  covered  panel  (V)  option  should  be  used  in  most  cases  at  the  wing-body 
intersection,  particularly  if  a  portion  of  the  wing  stub  is  defined  as  part  of 
the  body,  or  if  lateral  corners  with  ,  GO0  are  adjacent  to  the  covered 

panel. 

A  new  section  is  required  to  change  the  number  of  input  points  per  row  or 
to  introduce  an  inlet-type  step.  At  each  new  section,  the  spline  fit  starts 
anew,  and  the  slopes  and  ordinates  may  be  different  from  those  at  the  end  of 
the  previous  section. 

Normally,  the  slope  of  the  lateral  spline  fit  is  +  90°  at  the  first  and 
last  points  in  a  row.  If  this  is  not  desired,  these  points  must  have  a  CY 
point  type.  Also,  the  first  and  last  points  in  each  lateral  row  must  lie  in 
the  plane  of  symmetry  (  Y<10”®x  body  length).  All  longitudinal  lines 
defining  lateral  boundaries  of  panels  must  have  positive  length  (  i  X  >  2  x 
10'6X  body  length). 

The  fore  and  aft  panel  boundaries  (lateral  input  rows)  have  a  maximum 
slope  limitation  because  of  the  marching  scheme  used  to  solve  the  non-linear 
equations.  The  marching  scheme  requires  that  succeeding  rows  not  have  any 
effect  on  the  current  row.  Thus,  the  characteristics  from  the  back  of  a  row 
must  not  trace  out  of  the  back  of  the  row.  This  means  that  the  lateral  slope 
must  always  be  less  than  /f1t  £l/v/l  +  where  Mt  is  the  true 

local  llach  number  with  a  minimum  value  of  1.1  and  e  is  the  streamwise  surface 
slope.  To  meet  this  requirement,  all  lateral  lines  should  have  a  slope  less 
than  0.450/  J 1  + 


i 
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When  wing  data  are  input,  the  airfoils  can  be  defined  by  either  of  two 
formats  of  area  rule  type  input  (airfoil  ordinates)  or  by  streamwise  surface 
slope  input.  When  surface  slope  input  is  used,  care  must  be  taken  with  the 
sign  of  the  slope.  The  input  is  -1//2  dt/dx,  so  the  sign  is  negative  when  the 
thickness  is  increasing,  and  positive  when  the  thickness  is  decreasing. 

A  restart  capability  is  available  to  prevent  having  to  duplicate  long 
computations  in  cases  such  as  encountering  the  CP  time  limit,  a  program  bomb, 
or  to  evaluate  effects  of  goemetry  changes  near  the  aft  end  of  the  body.  To 
use  this  capability,  the  file  TAPES  must  be  saved  from  the  initial 
calculation.  Then  the  calculations  can  be  restarted  at  the  beginning  of  any 
row  for  which  all  preceding  row  calculations  were  completed.  E.g.,  if  the 
initial  calculation  was  completed  through  row  6,  the  case  can  be  restarted  at 
any  row  from  2  through  7.  The  Mach  number  must  be  the  same  ana  the  geometry 
through  the  preceeding  row  must  be  unchanged.  The  program  is  signaled  to  use 
the  restart  capability  by  reading  the  input  Mach  numbers  ano  finding  tne 
second  "Mach  number"  is  *  -2.  The  calculation  starts  at  the  'n'th  row  if  the 
second  "Mach  number"  is  (-n).  The  old  TAPE8  file  must  be  attached  to  the  job 
with  the  local  file  name  of  TAPES. 


7.0  ARRAY  L1M11S 


a) 

b) 

c) 

d) 

e) 

f) 


NPANL  • 

NROWS  * 

NPL  S 
N YD  I  SC  S 

MLDSC 

NSECS 


200,  where  NPANL  =  actual  numDer  of  panels  +  1  extra  idur.imy)  per 

row. 

20,  where  NROWS  =  numbers  of  rows  of  panels  (is  less  than  number 
of  input  stations  by  the  number  of  sections). 

19,  where  NPL  is  the  number  of  actual  panels  in  any  row. 

33-NPl,  where  NYUISC  is  the  number  of  Y,  or  lateral, 
discontinuities  in  the  row. 

s  80,  where  MLDSC  is  the  total  number  of  Y  (lateral)  slope 
discontinuities. 

s  b,  where  NSECS  is  the  number  of  sections. 


g)  Total  number  of  input  points  $  400. 
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8.0  OUTPUT  DATA  FROM  ANALYSIS  ROUTINE 


The  output  data  consists  of  the  following  categories  in  order: 

a.  Date  and  time  of  run 

b.  Title  cards 

c.  Body  geometry  input  -  literal  and  interpreted,  side  by  side 

d.  Wing  input  data-interpreted-if  used. 

e.  Reference  area  card,  Mach  numbers  -  interpreted 

f.  Summary  of  sections,  rows,  points 

g.  Wing-alone  results  -  if  applicable 

h.  Panel  geoemetry 

i .  Body  results 

j.  Results  for  wing  with  effect  of  body  -  if  applicable 

Only  items  (g)  through  (j)  above  will  be  discussed,  because  (a)  through  (f) 
are  self-explanatory. 

The  wing  alone  results  include  chordwise  data  at  each  span  station  and  the 
integrated  drag  coefficient  of  the  wing  based  on  the  specified  reference 
area.  The  span  stations  (y's)  are  relative  to  the  y  of  the  root  chord,  RUTY . 
At  every  5%  of  the  local  chord,  values  of  upper  surface  ordinate  (Z/C)  as  a 
fraction  of  local  chord,  local  slope  (ALP=  -dz/dx)  and  pressure  coefficient 
(CPU)  are  printed.  The  CDC  value  is  the  local  section  drag  coefficient  times 
the  local  chord,  divided  by  the  average  chord  of  the  exposed  wing. 

The  panel  geometry  print  gives  geometric  data  on  each  panel  at  the  four 
corners  and  at  the  average  x  on  each  side  of  the  panel.  The  RDD  is  a 
curvature  measure;  i.e.,  the  second  derivative  of  radius  with  respect  to  X. 

The  ordinates  X,  Y  and  Z,  the  derivatives  YU  (dy/dx)  and  ZD  (dz/dx),  the 
streamwise  SLOPE  and  the  outward  surface  normal  THETA  are  printed  at  each 
point  from  the  parametric  cubic  surface  patch  definitons.  The  analysis 
routine  uses  a  circular  arc  approximation  at  each  X  to  the  surface  patch 
defined  shape.  The  1/RADIUS  and  ARC-THETA  are  the  reciprocal  of  the  circular 
arc  radius  and  the  0,  respectively,  associated  with  the  circular  arc 
approximation.  Theta  is  the  direction  of  the  projection  into  an  X  -  constant 
plane  of  the  local  surface  outward  normal,  and  :  =  0  is  parallel  to  the 
positive  y  axis. 
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The  body  results  print  the  calculated  pressures  and  flow  conditions  at 
each  control  point.  Control  point  calculations  are  made  at  the  beginning 
(front)  of  the  first  row  and  for  other  rows  when  required  due  to 
discontinuities.  Control  point  calculations  are  always  made  at  the  middle  and 
back  of  each  row.  X,  Y  and  Z  are  the  ordinates  of  the  control  point.  Control 
points  are  located  20  percent  of  the  panel  width  away  from  the  side  edges, 
except  at  the  leading  edge  of  finite  width  panels  where  they  are  located  10 
percent  of  the  panel  width  away  from  the  edge.  THETA  is  the  lateral  outward 
normal  direction  as  defined  previously;  it  is  not  a  meridian  angle  except  in 
very  simple  cases.  EPS  is  the  local  strearnwise  slope.  CP  is  the  pressure 
coefficient.  M-LOCl  is  the  true  local  Mach  number.  PHIX,  PHIL  AND  PHIN  are 
the  perturbation  velocity  ratios  to  freesteam  velocity  in  the  X,  lateral  and 
normal  directions.  PHIL  is  positive  in  the  downward  or  increasing  panel 
number  direction.  V/VZ  is  the  ratio  of  the  total  perturbed  local  velocity  to 

freestream  velocity.  S3-EDGE  is  the  source  strength  at  the  adjacent  lateral 
edge  of  the  panel;  this  is  the  only  printed  value  which  does  not  correspond  to 
the  control  point  location.  BETA-C  is  the  correlated^  .  PH  I Y  and  PHIZ  are 
the  perturbation  velocity  components  in  the  Y  and  Z  directions.  Two 
additional  numbers  may  be  printed  at  the  right  side  of  the  page  without 
labels.  The  first  one,  which  will  include  a  decimal,  is  a  total  pressure 
ratio  which  arose  from  a  blunt  nose  or  corner  solution  and  is  attenuated  as 
the  solution  marches  aft.  The  second  number  is  an  integer  and  tells  how  many 
iterations  were  computed  for  the  point  or  group  of  points.  The  maximum  number 
of  iterations  allowed  is  five.  When  calculations  for  a  row  are  completed,  the 
contributions  of  each  panel  (one  side  only)  to  drag,  lift  and  pitching  moment 
are  printed  (DEL-CD,  DEL-CL,  DEL-CM).  Then  the  total  drag,  lift  and  moment 
coefficients  for  both  sides  of  the  configuration  and  summed  over  all  completed 
rows  are  printed  as  CD,  CL  and  CM. 

The  results  for  wing  with  the  effect  of  the  body  are  in  the  same  format  as 
the  wing-alone  results  described  previously,  and  are  printed  after  all  body 
results  are  listed. 


9.0  RESOURCE  REQUIREMENTS 

The  central  processor  time  per  Mach  number  presently  varies  from  1  to  15 
seconds  per  panel  on  a  CDC  6600.  If  the  input  geometry  is  axisymmetr ic,  the 
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routine  only  calculates  one  solution  per  streanwise  location  and  assigns  the 
resulting  pressures  to  all  lateral  panel  points,  thus  significantly  reducing 
the  computations  required.  Approximately  200,0003  words  of  memory  are 
required  for  the  analysis  routine  with  segmentation. 
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